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INTRODUCTION 


1  . 

An  infrared  (IR)  detection  system  determines  the 
presence  or  absence  of  a  target  in  a  particular  area  of 
space  by  processing  the  received  1R  energy  from  that  area. 
The  target  must  be  found  amid  ever-present  background 
radiance  and  any  internal  noise  generated  by  the  detection 
system.  In  this  thesis,  one  type  of  system,  called  a 
nutating  optical  system,  is  considered. 

Infrared  detection  is  a  combination  of  spatial  frequency 
processing  and  temporal  processing.  Spatial  frequency  process 
ing  is  a  means  of  achieving  a  modulated  output  from  the  IR 
radiance  passing  into  the  optical  system  so  that  information 
may  bo  easily  extracted  [Ref.  14].  To  extract  this  informa¬ 
tion,  some  form  of  temporal  processing  is  used.  In  a 
nutating  optical  system,  the  optical  axis  rotates  about  an 
axis  perpendicular  to  the  image  plane  in  such  a  way  that  a 
point  image  traces  out  a  circle  in  the  image  plane.  The 
spatial  frequency  filter  is  a  piece  of  IR  sensitive  material 
in  the  image  plane  which  produces  a  voltage  output  that  is 
amplified  and  fed  Anto  a  temporal  filter  for  processing. 

A  nutating  optical  system  is  frequently  used  in  a  missile 
seeker  because  it  produces  a  signal  when  there  is  no 
tracking  error  and  the  rotating  optics  can  be  made  into  a 
gyro  that  provides  the  inertial  reference. 

This  work  is  a  synthesis  and  extension  of  two  recently 
published  articles  in  the  literature.  The  first,  by 


Hargcr  [Rei‘.  3],  formulates  tlic  detection  theory  of  a 
known  signal  in  background  and  white  noise.  Since  in  1R 
target  detection  the  signal  is  never  known  completely, 
llargcr's  work  is  extended  to  cite  case  where  the  target 
amplitude  and  position  arc  unknown.  However,  a  known 
target  shape  is  assumed  us  well  us  a  priori  statistics  of 
amplitude  and  position.  Since  detection  is  most  difficult 
when  the  amplitude  is  small,  a  test  (called  the  threshold 
detector)  is  derived  that  is  optimum  for  weak  signals. 

To  specify  the  form  of  the  threshold  detector  and  to 
calculate  it:  per formance ,  an  eigenvalue  integral  equation 
must  be  solved.  Ordinarily,  this  is  a  difficult  undertaking 
but  the  solution  becomes  trivial  if  t lrc  covariance  function 
of  the  background  process  is  periodic.  This  is  exactly 
the  condition  realized  for  the  background  process  out  of 
a  nutating  detector  provided  that  the  phase,  which  is 
physically  unimportant,  is  averaged  out. 

Samuclsson ' s  results  [Ref.  11]  specify  the  form  of  the 
signal  and  background  coefficients  in  terms  of  the  system 
parameters  and  signal  and  background  radiance.  These 
derivations  were  first  published  in  a  Swedish  internal 
report  [Ref.  10],  then  re-dcrived  in  English  by  Yoshitani 
[Ref.  6],  The  latter's  derivations  arc  included  in  the 
appendices . 

Samuclsson's  equations  have  been  applied  to  a  rectangu¬ 
lar  IR  detector  assuming  a  Wiener  spectrum  for  the  random 
background  and  the  coefficients  have  been  calculated. 
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I  I.  SPAT!  A!.  PROCESSING _ 

Tlic  spatial  frequency  filtering  portion  of  the  detection 
system  is  comprised  of  a  lens  system  and  an  IR  sensitive 
detecto*  located  in  Tt  c  image  plane,  the  focal  plane 
of  the  optics.  To  adequately  derive  information 

from  the  IR  radiation,  the  lens  system  and  detector  must 
be  of  a  nature  that  causes  the  voltage  output  to  be  modulated 
in  a  way  that  information  about  a  target  within  the  field 
of  view  can  be  processed.  In  a  practical  nutating  optical 
system,  the  lens  system  moves  and  the  detector  is  stationary. 
However,  for  conceptual  purposes,  this  is  equivalent  to  a 
stationary  lens  systen  and  nutating  detector. 


A,  TARGhT  SPECTRUM 

The  object  to  be  detected,  the  target,  is  located 
jn  the  object  plane  a  distance  R  from  the  lens  and  parallel 
to  the  image  plane  as  shown  in  Figure  1.1.  The  angle  fro"-' 
the  perpendicular  axis  connecting  the  planes  to  the  target 
coordinates  measured  in  the  x -direction  is  called  x  and  is 
measured  in  milliradi ans ;  the  angle  in  the  y-diiection  is 
y.  The  target  is  considered  to  be  close  to  the  optical 
axis  and  at  a  sufficient  distance,  R,  such  that  small  angle 
approximations  may  be  made, 


x  5  tan  x 
'/  "  t£,n  '/ 


(2.1) 
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image  plane 


Consider  a  target  with  a  radiance  distribution 
2 

N  (r)  (W/cm  -sr)  which  is  the  radiant  power  into  a  unit  solid 
angle  per  unit  projected  area  of  the  sources.  r  =r(x,y)  is 
a  vector  whose  origin  is  the  perpendicular  axis  connecting 
the  coordinates. 

Assuming  no  atmospheric  attenuation  (this  can  be 
included  in  N(r)  if  it  exists)  and  perfect  optics,  the  power 
distribution  in  the  image  plane  is 

P i ( r )  -  N(r)  *  \Q  (w/sr)  (2.2) 

where  Aq  is  tnc  effective  entrance  area  of  the  optics. 

Dividing  by  Aq ,  one  can  see  that  the  radiance 
distribution  in  the  image  plane  is  numerically  equal  to  the 
radiance  distribution  in  the  object  plane. 

However,  because  the  optics  is  not  perfect,  a 
point  object  given  by 

Hp  6 (r)  -  Hp  5 (x)  5  (y)  (w)  ,  (2.3) 

2 

(where  Hp  (l\'/cm“)  is  the  irradiance  at  the  optics  received 
from  the  point  source)  is  imaged  as  a  power  distribution 

Hp  •  fQ(r)  (K/cm2  Sr)  (2.4) 

where  fg(r)  (sr  2)  is  the  point  spread  function  of  the 
optics.  Irradiance  is  defined  as  the  radiant  flux  indicent 
on  a  surface  of  unit  area  [Ref.  5]. 

Therefore,  the  radiance  distribution,  N(r),  in  the 
object  plane  will  be  imaged  as 

N'(r)  =  Jn(t-s)  fQ(s)  d2s  (W/cm2Sr)  (2.S) 
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2 

where  d  s  *=  dxdy  and  the  integration  is  over  the  entire 
image  plane. 

The  detector  is  a  space  filter  in  the  image  plane 
which  may  be  a  moving  retical  or  simply  an  aperture  across 
which  the  power  is  integrated. 

The  power  in  the  image  plane  incident  on  the  detec¬ 
tor  is 

H-  j N'(r)  t(r)  d2r  (K/cm2)  (2.6) 

where  r  (r)  -  t(x,y)  is  tlic  transmittance  of  the  detector. 

When  the  dctcctoi  moves  as  a  function  of  time,  the 
power  incident  on  the  detector  becomes  ;>  function  of  time, 

HCt)  -  j N '  (r)  r (r - p (t) )  d2r  (W/cm2)  (2.7) 

where  p(t)  describes  the  movement  of  the  detector  in  the 
image  plane. 

A  nutating  detector  moves  circularly  in  the  image 
piano  but  its  orientation  remains  fixed,  i.e.,  the  coordin  te 
system  (x.’  ,y')  shown  in  Fig.  1.2  does  not  rotate.  The  nut 
tion  rad. us  measured  with  respect  to  the  optical  axis  is 
given  by  p. 

Assuming  a  nutating  optics  and  no  relative  motion 
between  the  target  and  the  optics,  the  radiance  on  the 
detector  will  be  periodic  and  thus  can  be  expanded  in  a 
Fourier  series 

«  j  mo  t 

H(t)  *=  l  Hn  c  0  0  t  <  T  (2.8) 

u 
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Replacing  H(t)  by  equation  (2.7], 


i  fT  -jniijt 

Hn  T  J  dt  e  d  J  N’(r)  (r-p(t))  (2.9) 

The  motion  describing  function,  p(t),  for  circular 
nutation  is 

p(t)Ep(cosw0tlsinoJ0t),  (2.10) 

where  =  2  rr/T  (rad/sec)  is  the  radian  frequency  of  nutation 
and  T  is  the  period  for  one  nutation. 

Substituting  (2.10)  in  (2.9),  II  can  be  shown  (see 
Appendix  A)  to  be 

Hn  «=  jn  Jn  '  (k)  T*(k)  e'jn*  2n(2.V^  d2k  (2  11) 


Figure  1.2.  Circular  Nutation  in  the  Image  Plane. 


where  k 


=  (k  ,  k  ),  two-dimensional  spatial  Courier 

-  x  y 

transform  coordinates  [Ref.  8], 

N'(k)  ■  transform  of  radiance  distribution  in  the 
image  plane, 

i*(k)  ■  complex  conjugate  of  the  transform  of  the 

optical  transmittance  function  in  the  detector 
coordinate  system, 

❖  “  tan'1  ky / kx , 

p  =  nutation  radius, 


and  J  (2)  is  the  nth  order  Bessel  function  of  the  first 
kind.  It  has  been  assumed  that  the  image  plane  is  very 
large  compared  to  the  radiance  distribution  from  a  target 
so  that  the  image  plane  may  be  considered  infinite  for 
mathematical  purposes. 


B.  CORRELATION  OF  BACKGROUND  NCISE 

The  background  power  incident  on  the  detector  is  from 
a  sample  background  scene  where  a  scene  is  a  two-dimensional 
random  process  characterized  by  the  radius  vector  r  =  (x,y) 

from  the  optical  axis. 

Let  N(r)  be  a  sample  scene  radiance  distribution  on  the 

object  plane  from  an  ensemble  of  scenes  which  lias  been  given 

suitable  probability  structure.  The  background  power  incidcn 
on  the  detector  is 

B(t)  -  /  N"  (r)  x  C  r  -  p  C  t )  )  dZr  .  (2.12) 

The  covariance  function  is  defined  by  taking  the  expected 


value 


E{ [B(t)>  -  E{fl(t)J] (B(u)  -  L{H(u)}]  -Rfi(t,u)  (2.13) 

where  E  denotes  the  expectation  over  the  ensemble  of  scenes. 
Because  the  detector  is  nutating,  the  random  process,  B(t), 
will  not  be  stationery.  However,  if  the  detector  and 
background  do  net  move  with  respect  to  each  other,  then 
Kg  ( 1 1 » 1 2 )  is  doubl>'  periodic: 

RgCtj  ♦  MT,  t2  +  NT)  -  RBCt1,t2)  (2.14) 

where  M  ai.d  N  are  integers. 

It  ran  be  shown  [Ref.  7]  that  the  form  of  the  covari¬ 
ance  of  a  doubly  period  process  is 


RgCtj.tO 


I 

m,n 


JmVi 


mn 


c‘ jnuot2  . 


(2. IS) 


Because  the  nutation  is  periodic,  over  the  ensemble  of 
scenes,  there  is  little  significance  in  the  starting  time  t. 
Therefore,  considering  t  as  a  random  variable  that  is 
uniform  over  the  nutation  period,  one  can  "average  it  out" 
of  t!  c  covariance  function  as 


Rb(i)  -  \  /  RB(t+T,t)  dt  . 


(2.16) 


Substituting  (2.16)  in  (Z.15), 

1  -T  j  (m-n)vJ0t 


*8^  *  1  6mn  ' 

m  ,n 


w  J  C 
1  0 


dt. 


(2.17) 
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Using  the  relation 


i  T  J(u-v)unt 

6uv  ■>  ijr  /  e  at  ,  (2.181 

1  0 

the  correlation  function  becomes 

®  ;inu  r 

RbCt)  -  r  6n  c  u  .  (2.193 

n 

Equation  (2.19)  shows  that  B(t)  can  be  considered  a 
wide  sense  stationary  process  and  that  it  is  periodic  in 
the  mean  square  sense;  i.e., 

Rp ( t )  “  Rb(T+t)  for  any  t  (2.203 

The  power  spectral  density  is  the  temporal  Fourier 
transform  of  Rb(t).  Taking  the  transform  one  has 

00 

Sb(«3  1  2*  Z  -Uu-no.0).  (2.213 

n 

The  coefficients,  3  ,  can  be  related  to  the  optical 
system  parameters  by 

Bn  “  /iT  I  2  l  ^'o  Ci2)  1  2  VIP  Jn2  (ZlTp/\  ^  (2‘22) 

where  ( k 3  is  the  transform  of  the  point  spread  function, 
and  Wg(k3  is  the  Wiener  spectrum,  the  transform  of  the  back 
ground  correlation  function.  This  derivation  was  first  mad 
by  Samuclsson  [Ref.  10]  in  an  internal  Swedish  report  then 
again  by  Yoshitani  [Ref.  &].  Yosluinti’s  derivation  has 
been  included  as  Appendix  B. 


liquation  (2.  ID)  implies  that  Hie  background,  B(t),  can 
be  represented  by  a  Fourier  series  with  uncorrelatcd 
coefficients  [Ref.  7].  Thus, 

"  j»unt 

B  (t)  c  I  b  c  U  (2.23) 

n  n 

where  the  convergence  is  in  the  mean  square  sense. 

The  coefficients,  b  ,  satisfy 


BCbn) 


E(il(t))  n 


n  f  0 


(2.24) 


Moreover,  they  arc  uncorrelatcd 


C.  WHITE  NOISE  COMPONENT 

In  an  actual  system  there  will  be  a  certain  amount  of 
internal  noise  generated  in  the  detector  and  preamplifier 
which  is  additive  to  any  background  noise  or  signal.  In 
the  frequencies  of  interest  the  primary  source  of  this 
noise  is  thermal  agitation.  This  type  of  noise,  called 
Johnson  noise,  is  assumed  to  be  white  and  Gaussian. 


Let  W(t)  be  Gaussian  white  noise  with  power  spectral 
N  7 

density  ( \r  *"  / 1 1  )  .  The  power  spectrum  at  the  output  of  the 

preamp  with  bandwidth  U  is  shewn  in  Figure  1.3. 


Mag  . 
(V/HO 


III.  TEMPORAL  PROCESSING 


The  output  voltage  from  the  detector  is  temporally  pro¬ 
cessed  to  determine  the  presence  of  a  target.  The  type  of 
processor  is  determined  by  applying  statistical  detection 
theory  to  The  known  beho\  lor  of  the  nutating  system  and  var¬ 
ious  realistic  assumptions  about  the  system  and  the  target. 

An  optimum  processor  in  the  sense  of  small  signal  amplitudes, 
called  a  threshold  detector,  is  developed  in  this  section. 

Harger  [Ref.  3]  derived  the  detection  theory  for  the 
known  signal  ease  in  the  presence  of  additive  background  and 
white  noise  components.  In  IR  detection,  the  target  ampli¬ 
tude  and  position  are  usually  not  known  completely.  The 
threshold  detector  is  an  extension  of  Hanger’s  v. ork.  which 
accounts  for  the  unknown  amplitude  and  position. 

The  form  of  the  optimal  detector  under  the  Neyman- Pearson 
criterion  will  be  derived  for  multiple  observations.  It  i; 
assumed  that  the  image  plane  and  object  plane  remain  parallel  .  \d 
stationary  with  respect  to  each  other  during  the  observati  ns. 

The  detection  of  the  target  under  the  Neyman-Pearson 
criterion  becomes  a  problem  of  hypothesis  testing.  The 
task  of  the  detector  is  one  of  choosing  between  two 
hypotheses,  Hq  that  only  noise  is  present  and  that  in 

addition  to  the  noise  there  is  a  target  present.  The  design 
of  the  detector  is  one  that  permits  the  correct  choice  of 
hypothesis  (a  detection)  with  maximum  probability  of 
detection,  Qd,  in  a  fixed  probability,  Q^,  of  choosing 


when  Hq  is  true  (a  false  alarm).  The  structure  of  the 
Ncyman-i'curson  criterion  requires  forming  a  likelihood  ratio, 
which  is  compared  to  a  known  threshold,  V  ,  specified 
for  a  given  false  alarm  probability.  If  the  likelihood 
ratio  exceeds  the  threshold,  hypothesis  is  assumed  true 
and  if  not,  hypothesis  Hq  is  assumed  true.  The  logic  process 
may  be  written 


10 


(3.1) 


where  the  threshold,  V  ,  is  erived  from  the  expression 

co 

Qj.'A  e  /  P  (Zt  Hx)  dZ  (3.2) 

Vt 

and  Z  is  a  sufficient  statistic  of  the  received  data. 


A.  DETECTION  OK  A  KNOWN  SIGNAL 

The  data  for  processing  is  a  set  of  M  functions  of  time 
(Zn(t),  where  each  function  is  the  output  of  the 

preamp  during  one  nutation  (0,T).  T  is  the  nutation  period. 
The  output  may  be  represented  as 

ZnCO  -  Sn(t)  +  B(t)  +Nn(t);  0  <_  t  <_  T ,  n-1 . M  (3.3) 

S  (t)  represents  the  output  due  to  a  target.  The  component, 
B(t),  represents  an  additive  background  and  is  a  sample 
function  of  a  random  field  (B) .  The  relative  si2c  of  the 
target  in  the  field  of  view  is  assumed  sufficiently  small 
that  the  background  is  additive  while  B(t)  is  assumed  to  be 
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the  sample  function  during  each  of  the  M  nutations  observed 

Nn(t)  represents  the  additive  noise  component  due  to  the 

internal  noise  of  the  detector-preamp  combination.  Tt  is  a 

sample  function  of  a  white  Gaussian  random  process  ( W ^ )  of 

2 

zero  mean  and  power  spectral  density  N/2(V 

The  two  hypotheses  for  the  likelihood  ratio  may  be 
written  as 

Hy  <Zn(t  )  -  Sn(t  )  +  B(t  )  +  Nn(t  )  ;  n-l,...,M} 


Hq  :  { 2n ( t  )  -  B(t  )  ♦  Nn(t  )  ;  ne 1 . M 

Denote  the  likelihood  ratio  as  A,n({Z  })  where  the  like 

1 U  n 


lihood  ratio  is  defined  as 


Ai0CUnn 


P(tZn>iH,) 


'n ;  I  0  ■ 


C3.5) 


The  likelihood  ratio  is  complicated  bv  the  fact  Z  and 

n 

Zm  are  not  independent  because  of  a  common  noise  component 
B(t).  To  simplify  the  derivation  Harger  introduced  an 


auxiliary  hypothesis 


H2:  <?n  '  Nn  ;  n‘l . M)- 


C5.6) 


Then  the  likelihood  ratio,  A ,  can  be  computed  using 
the  chain  rule  for  likelihood  ratios: 

paying  pc{zn)|Hi)/P({zn}|H2) 


'n  »  0 


'n  I  0' 


P  u  y  | 


A 1 2  ^  A  0  2 
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If  the  sample  functions  for  the  background,  B ,  were  non- 


random,  Z.(t)  und  Zj  (t) ,  i  i  j,  would  be  mutually  indepen¬ 
dent  because  of  white  noise.  Accordingly,  if  B  were  fixed, 
one  can  write  the  conditional  likelihood  for  the  n1^ 
observation  C  { 2  n  }  ]  B ")  .  Moreover,  the  likelihood  ratio 

for  M  observations  would  be  the  product  of  the  likelihood 
ratios  for  each  single  observation, 

nc  1 


The  conditioning  is  then  removed  by  averaging  over  the 
ensemble  of  backgrounds 


A  i  -> 

X  <- 


(Ur>) 

-  1 1 


Eb|a„((2,HB)  1  . 

A4.  “*  44  —  J 


Since  the  noise  is  white  Gaussian,  the  conditional 

likelihood  ratio  is  that  for  detecting  a  known  signal 

S  +  B  in  white  Gaussian  noise  and  has  the  well  known  form 
n  n 

([Ref.  IS] ,  page  253) 


(fzn>|B) 


exp 


1 

N 


M 

E 

n*- 1 


T 

f 

0 


1' 

/  2_(t)(S  n(t)  *  B  ( t) )  dt 
0  n 

?  ] 

(Sn(t)  +  BCt))-  dt 


which  may  be  rewritten  as 


A12  Zn} I c  exp 


M 


T 


jf  S  '  Vt}  dt 

n=l  0  n  n 


M 

E 

n=  1 


S7  (t) 
n  v  ‘ 


dt 


exp  u  (b) 


(3.14) 


23 


a(b) 


(3.15) 


"  K  z  *  (zn(t^ 
N  n=J  0  n 

To  line!  Aio({Zr}j  , 


+  Sit))  B(t)dt  -  S-  /  B  1 1 )  d  t . 
n  N  0 


one  must  evaluate  Eg  ( A11^ ^  ( {  }  |  B)  ] 


where 


M  T 

EB[Ani2({Zn}  !PJ^  =  exP  1  f  Zn^  Sn(t)  dt 

n=l  0 

(3.16) 

M  T 

-  V  Z  f  S  ‘(t)ilt)  •  E,.  {exp  a  (b )  }  • 
n=l  0 


Assunr'ng  the  background  is  Gaussian  the  expectation  can  be 
evaluated  by  r  presenting  the  background  in  a  Karhuncn- 
Locvc  (K-L)  expansion 


B(t) 


Z  b, 


$k(t)  .  0  <  t  <  1 


(3.17) 


where  $k(t)  is  the  eigenf  notion  cor i  esponding  to  eigenvalue 
of  the  homogeneous  integra1  equation 

T 

^k9k(t)  /  Rg(t,u)  <f>k(u)du  ,  0  <  t  <  T.  (3.18) 

BgCtjU)  is  the  covariance  function  of  background 

RB(t>u)  l:{[B(t)  -  mB(t)][B(u)  -  nR(u)]}  (3.19) 

where  mR(t)  is  the  mean  of  the  background. 
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Then  Eg(expa(b))  can  be  calculated  in  a  straightforward 
manner.  After  thi  lengthy  calculation  shown  in  Appendix  C 


In  En  (exp  a(b])  =  -  l  E  In  (1  +  2  X,  ) 
B  l  k  X 


-  T  M  T  M 

♦  w  /  dt  Z  (Z  (t)  -  S  (t))  f  du  E  (Zm(u)-Sm(u)3q(t,u) 
0  n=  1  0  m=l 


.  T  -  M 

1  /  dt  {mR (t)  •  £  E  (Z  (t) -S  (t) ) }  h(t)  (3,203 

2  0  B  n=l  n  .  n 


where  q(t,u3  satisfies  the  integral  equation 


q  ( t ,  u  3  +  /  RB(t,i)q(t  ,u3rlT  =  RjjCt.u)  0  <  U,  t  <  T  (3.213 
and  h(t)  satisiies 


h  (t)  +  f  Rg  ( t ,  1 3 h  (x 3  dt  =  mH  (t)  0  <  l  <  T, 


0 


(3.223 


q(t,u3  is  the  impulse  response  of  the  filter  that  gives 
the  minimum  mean  square  error  estimate  E  of  the  backgrouni 
when  the  mean  is  zero  [Ref.  15].  The  additive  white  noise 
is  reduced  in  amplitude  by  1/M.  Using  the  K-L  representa¬ 
tion,  one  can  also  write 


q(t,u)  = 


N  > 

TET  Ak 

N 

m  +  h 


«k(t)  <>k(u)  0  i  u»  t  <  T  (3.233 


a  i 


Using  equations  (2.93,  (3.103,  (3.14}  and  (3.20)  and 
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-i*** 


**  iSJ 


observing  that 


A02((Z„»  •  A12tl2n»lSn  ■  0 


(5.2-1) 


it  can  be  shown  (sec  Appendix  C)  that  the  likelihood  function 
for  hypotheses  and  Hq  is 


A10«Zn» 


M  T 

exp  {  I  /  Z  (t)q_(t)dt  +  K' } 
n-1  0  n  n 


(3.25) 


where 


qn(t) 


N  [Sn^^ 


1  M  T 


T*  M 

i  /  I  S  (u) q(t ,u) du] 

w  o  £=1 


TT  M  .  M 


(3.25a) 


FT  £  /  C  t)  dx:  +  //  E  t  S(t)S,  (u)q(t,u)dtdu 

n=l  0  n  ,NM  00n»l  £-1  n  * 


7  M  T 

-  |  l  /  S  (t)h(t)dt 
m  n=  1  0 


(3.25b) 


The  logarithm  is  a  mcnotonic  function.  Thus  it  offers 
no  loss  of  information  but  greatly  simplifies  the  calcula¬ 


tions.  In  the  remaining  work 


A'l0  ~  ln  A10  '{Zn^  4 


(3.26) 


Rewriting  equation  (3.25), 


A'lO(Un})  =  L  ([Zn})  +  K’ 


(3.27) 


where 


U 
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M  T 

L(  Z  )  “  l  f  Z  (t)q  (t)dt  . 
n  n-1  0  n  n 


K'  does  not  depend  on  the  input  data  so  it  nay  be  included 
in  the  threshold.  The  hypotheses  test  now  becomes 


L (  Z  )  „  Threshold  .  (3 

n  < 

H0 

It  is  also  possible,  as  Harger  lias  shown,  to  derive 
tests  similar  to  (3.28)  without  making  the  Gaussian  back¬ 
ground  assumption.  In  particular,  if  the  background  fluctu 
ates  many  times  over  (0,  T)  such  that  a(h)  o 1  (3.15)  is 
approximately  Gaussian,  then  the  test  is 


M  T 


l  /  Z  (t)g  (t)dt  Threshold 

_  i  n  n  n 


nc  1  0 


where 


gn(t) 


T  M 

s  1S„U)  -  “  '  13 


Since  by  Mercer's  theorem 


Rg(t,u)  -  I  ^kCk(t)  <p£(u) 


0  <  t,  u  <  T 


on  comparing  q  ,  and  (3.21),  it  is  seen  that  the  above 
model  approaches  the  Gaussian  model  if 


7R  >;>  m£x 


The  above  condition,  which  implies  a  weak  background,  is 
sufficient  but  by  no  means  necessary.  The  two  models 
approximate  each  other  whenever  the  second  term  of  and 
gn  are  nearly  identical.  The  point  is  that  while  the 
Gaussian  background  assumption  may  not  fit  the  reality, 
a  processor  based  on  the  assumption  is  probably  not  far 
from  being  the  optimum.  Also,  the  Gaussian  assumptioT  allows 
one  to  compute  the  processor  performance. 


B.  DETECTION  OF  SIGNAL  WITH  UNKNOWN  PARAMETERS 

The  likelihood  ratio  derived  in  section  A  is  applicable 
only  if  the  signal  is  known  completely.  Such  is  seldom  the 
case  in  IR  target  detection  and  in  general  the  signal  is  a 
function  of  some  unknown  parameters.  Typically,  the  ampli¬ 
tude,  a, of  the  target  and  its  position  ro  in  the  image  plane 
are  the  unknown  parameters,  and  the  form  of  the  signal  can 
be  rewritten  to  include  these  parameters: 


S 


n 


Sn(t;a, 


Furthermore,  by  assuming  that  only  point  targets  are  of 
interest  and  including  the  assumption  of  fixed  image  and 
object  planes  over  the  number  of  nutations  of  interest, 
the  shape  of  the  waveform  will  be  the  same  during  each 
nutation.  However,  the  amplitude  is  allowed  to  vary  between 


successive  nutations. 


The  signal  waveform  may  now  be  written  as 


Sn(t;a,^o^  "  an  fCt.V  (3.31) 

where  an  is  the  signal  amplitude  and  £(t,ro)  is  the  signal 
shape  for  a  target  at  r  . 

Rewriting  hypothesis  under  these  assumptions 
gives 

Hl':  "  anfft,^o^  +  BCr)  +  NnCO;  n,cl . M 

0  <  t  <  T}  (5.32) 

The  likelihood  ratio  A 1Q  becomes  the  conditional 
likelihood  ratio  (( Zn> | a,rQ)  obtained  by  replacing 
S^(t)  by  anf(t,rQ)  everywhere  in  (3.25). 

Assuming  a  priori  statistical  knowledge  of  a  and  Tq  , 
and  their  independence,  the  likelihood  ratio  may  be 
averaged  over  the  respective  density  functions  to  produce 
an  unconditional  likelihood  ratio 

A10({Zn})  c  //•W{V^Vp(!)?CVd?dV  (3’55) 

Instead  of  averaging  as  above,  one  could  estimate  a  and  rQ 
by  the  maximum  likelihood  procedure  and  substitute  the 
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estimated  values  a  and  in  the  conditional  likelihood 
ratio  and  try  to  maximize  the  generalized  likelihood  ratio 

A 

Aio  ({Z^} | a, Tq) .  However,  this  procedure  is  not  followed 

here  because  the  maximum  likelihood  estimates  are  difficult 
to  obtain. 

Since  the  main  interest  is  in  detection  0f  the  weak 
target,  the  small  signal  case  is  of  primary  interest. 
Expanding  A,q ( {Z^ } j a , t  )  in  powers  of  signal  amplitudes 


M 


A10({Zn}|a'rn)  ‘  1  +  E.  an  7T  Ain((Z„}|a,rn)  I. 


n-1  n  10  n  a*=0 


(3.34) 

*  *n-l  W  Vt5i^,10((zn1l!.f|,)l:.0 


M  M 


Integrating  (3.34)  as  in  (3.33)  term  by  term  and 
retaining  only  the  first  non-zero  term  that  depends  on  the 
data  provides  a  statistic  which  is  optimum  for  weak  signals 
and  is  called  a  threshold  detector  [Ref.  4], 

The  evaluation  is  made  easier  by  the  fact 


Aiot(V  Ij-V  'a70 


517  ln  AJO(t:n)|?-ro)ia.0 


which  results  from 

A({zn}|o,T0)  -  1  . 


(3.35) 
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The  first  non-zero  term  becomes 

Tjp  ln  A10"VLa’;0>laM)  ■  H  ^nCt)f(tJ0)dt 


T  T 

l  f  2  ft) dt  /  qCt,u)f (u,r0)du 
N  0  n  0 


T 


-  /  f(t,r  )h(t)dt 
0  u 


(3. Jo) 


Performing  the  integration  over  r Q , 

j  M  T  .  T 

4ln({zn  S  1  +  £  Z  a  /  Z  (t)  [f(t)  -  ;  q(t,u)  f(u)du] 
J.U  n  r*  n*i  n  o  n  0 


T. 

-/  fCt)  h(t)  dt 

o 


(3.37) 


where  the  bar  denotes  the  averaged  quantities.  The  first 
and  last  terms  are  independent  of  the  data  and  may  be 
included  in  the  threshold  value.  Moreover,  if  the  ampli¬ 
tudes  are  identically  distributed,  as  is  generally  the  case, 
the  threshold  detector  performs  the  test 


W 


M 

I 

n*=l 


H, 


;  Z  (t)  q(t)  dt 
0  n 


(3.38) 


where  W  is  the  threshold  and 

T 

q(t)  -  |  (f(D  -  /  q(t,u)  f(u)  du]  .  (3.39) 

0 
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Thus  the  threshold  detector  can  be  mechanized  either  as  a 
correlator  or  matched  filter  followed  by  an  accumulator. 

The  filter  function  q(t)  may  be  represented  in  a  series 
oxpansion  in  terms  of  the  eigenfunctions  of  Rg(t,u)  and  the 
coefficients  of  the  signal  expansion.  Expanding  f(t)  in 
terms  of  the  eigenfunctions: 


00 

fCt)  -  z  fk*k(t)  . 

Substituting  (3.40)  in  (3.39)  and  using  (3.23) 


q(f) 


♦  R(t)  • 


(3.40) 


(3.41) 


C.  EVALUATION  OF  PERFORMANCE 

The  Gaussian  assumption  was  made  in  deriving  the 
threshold  statistic;  therefore,  the  threshold  statistic 
itself  is  Gaussian.  This  means  that  the  probability  of 
detection,  Qd,  and  probability  of  false  alarm,  Q^,  are 
uniquely  determined  by  second-order  statistics,  e.g., 
the  means  and  variances  of  K  under  hypotheses  Hq  and  H^. 

Recall  that  the  threshold  statistic  is 

M  T 

W  -  I  /  Z  (t)  q (t )  dt.  (3.42) 

n-1  0  n 

The  means  and  variances  of  W  are: 
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M  T 


E(W|HQ)  -  £  /  q(t)  E(B(t)  ♦NJ1(t))-dt 


n-1  0 


M  T. 

£  /  q(t)  m  (t )  dt 


n-1  0 


B 


M  T 

E(W|H.)  -  £  /  q(t)  E[a  f(t,rfl)  +  B(t)  +  N  ft) ]  dt 

1  n-1  0 


M  T 

E(W|Hn)  *  £  /  f(t)  q(t)  dt 

0  n-1  0 


VAR  (W|Hq)  -  VAR  (Wjlip  -  VAR 


T  T  T 

NM  /  q  2  ( t)  dt  *  M2  /  q(t)  f  qCu)  RB(t,u)  dudt 


0 


0 


(3.43) 


Letting  pw(*lH^)  be  the  Gaussian  density  function  of 
W  under  which  is  specified  by  the  equations  above,  then 


«FA  '  (.  Pw'xll,0,dx 


(3.44) 


and 


<2d  *  MxlHPdx 

Wt 


(3.45) 


where  W  is  the  threshold  voltage. 

Equations  (3.43)  and  (3.44)  can  be  solved  to  yield 


(5.46) 


«J;A  ■  *cCx) 

and 

Qd  -  ♦c(x*d)  (3.47) 

where 

W  -  EfW|H  ) 

x  .  - -  , 

✓OT 

d  is  the  equivalent  s ignal - to -no i s c  ratio 

E(W|H  ] 
d  -  - - i_ 

/TO 


M  T. 

I  a  f  f(t)  q(t)  dt 
.  n-1  n  0 _ 


q2(t) 


dt 


2 

tr  ; 


Z - T - 

q(t)  /  q (u)  Ra ( t ,u)  dudt 
0  D 


(3.48) 


and  4>c(g)  is  the  complementary  error  function  of  the  kind 

1  *  -v2/2 

$,(g)  *■  -  /  e  dv  (3.49) 

d  vTF  g 

By  setting  a  desired  probability  of  false  alarm,  x 
may  be  determined  from  equation  (3.46).  The  probability 
of  detection  is  then  determined  by  solving  for  d  and  using 
equation  (3.47). 

By  substituting  (3.41)  and  (3.42),  and  using  (3.23), 
the  numerator  of  d  is 
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(3.50a) 


I  a  /  f(t)q(t)dt 
n-1  n  0 


l 

n*- 1 


X 


k 


while  the  denominator  of  d  (see  Appendix  D)  is 


ir 

7CT 


4 


Therefore  the  effective  signal-to-noise  ratio  is 


d 


1 


4- 


TFT 


(3.50b) 


(3.51) 


Knowledge  of  d  along  with  (3.46)  and  (3. ^7)  completely 
specifies  the  performance  of  the  threshold  detector. 

The  equivalent  signal-to-noise  ratio  for  a  target  whose 
amplitude  and  position  arc  known  is  easily  shown  to  be 


i  t 
M  nil 


i  [2 W 
n  yJT 


z 

k 


,(?0)|2 


1 

TFT 


1  *  IT  Xk 


(3.  2) 


By  letting  the  background  be  zero  it  will  now  be  shown 
that  d'  reduces  to  the  equivalent  signal-to-noise  ratio  for 
a  known  signal  in  white  Gaussian  noise  which  is  well  known 
(Ref.  4] . 

The  energy  dissipated  during  one  observation  in  a  resis¬ 
tor  of  1  ohm  if  the  signal  a  f(t,rQ)  is  the  voltage  across 
that  resistor  is 


(3.53) 


■  *n  jj  lfk'V>|J  • 

The  voltage  resulting  from  M  observations  then  is 

—  i  t  M  r* 

/E  ’  R  *.  /En  H  E  an  £  MV"  ■  (J.S4) 

n- 1  n*1!  k 

J 

Substituting  (3.54)  in  (3.52)  with  A  ^  ■  0, 

d'  -  /fT  fnT  -  I2E~  /M~  , 

\TT  V  N 

which  shows  the  familiar  result  that  detectability  depends 
only  upon  the  signal  energy  and  the  sagnai-to-noise  ratio 
increases  as  the  square  root  of  the  number  of  observations. 
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IV.  THRESHOLD  PFTHfTOR  FOR  NUTATING  SYSTEM 

In  an  infrared  detection  system  the  output  voltage  to 
be  processed  is  obtained  from  a  detector -prcampl if ier  com¬ 
bination.  The  irradiancc  on  the  detector  in  W/cm2  is  con¬ 
verted  to  volts  with  a  linear  scale  factor.  The  voltage 
output  can  be  written  as 

v(t,r0)  -  KH  (t  ,r0J  C 4  . 1 ) 

where  H(t,rg)  is  the  irradiance  on  the  detector  from  a  target 
at  position  rQ,  and  K  is  a  scale  factor  in  (V/W/cm2)  .  Also 
the  covariance  function  may  he  written  as 

RvRCt)  =  ElKB(tM)  KB (t) )  »  K2RB(rl,  (4.2) 

where  R(t)  is  the  irradiance  on  the  detector  due  to  a  partic- 
ular  background  scene.  The  output  voltage  due  to  this  scene 
is 

v g  «  KB ( t)  (4.3) 

The  problem  now  is  to  relate  the  detector -preamp  1 i fie v 
output  to  the  equations  derived  for  the  threshold  detector. 

Jt  was  shown  in  Part  III  that  an  optimum  processor  for 
weak  signals  was  a  threshold  detector  that  performed  a 
certain  linear  operation  on  the  data.  The  linear  filter  was 
specified  in  terms  of  the  eigenvalues  and  c-igenfunction^  of 
the  integral  equation  (3. IB),  where  the  kcrnal  R^ft.u)  is 
the  covariance  function  for  the  background. 
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In  Pert  '/I,  the  rovariancc  function  was  noted  to  be 
periodic.  Rewriting  (4.2)  in  terms  of  the  output  voltage 


7  ~  jVu  (t-u) 
RvB('.u)  ■  K  l  Sk  e 


(4.4) 


Substituting  (4.4)  into  integral  equation  (3.18),  the  solu¬ 
tion  for  the  eigenvalues  and  eigenfunctions  is  trivial  with 


JkV 


2 

v‘o  T 

k  v 


(4.5) 


where  T  is  the  nutation  period. 


1  he  signal  coefficients  weie  also  noted  to  be  periodic 
and  could  be  expanded  in  a  Fourier  series 


°°  j  k  uj  t 

v(t)  =  K  l  H,.(r  )e 
k  K 


(4.  ) 


Recognizing  the  similarity  hetv>een  equations  (4.6)  and  (3.  1)  , 

the  signal  coefficients  for  a  known  signal  which  is  expanded 


using  the  basis  set  above  become 


W  -  Kllk<V''T 


(4.7) 


For  the  threshold  detector,  the  signal  is  averaged  with 
respect  to  an  a  priori  density  function  for  r^.  Averaging 
can  be  done  with  the  coefficients  to  give 
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*k  "  1  p(^fk(,Vd'o 

*  K/T  /  1'j<(r0)p(^0)d^0  (4-8) 

-  K/T  Fk  . 

Using  equation  (2.10),  the  coefficients  Fk  can  he  written 

as 

\\  =  jn  /d^0P(r0)  /N'  (k,ro)x*(h)c^n^  Jp  [2^0  yfk^y  \  ]d2k  (4.9) 

The  transform  of  the  radiance  distribution  from  a  target, 
N'(k,rc)  is  a  function  of  the  target  radiance,  R(k),  the 
position  of  the  target,  r  ,  and  the  point  spread  function  of 
the  optics,  b  (k)  .  Rewriting  (4.9), 

Hk  -  jn  /d2kR(k)  t‘(k)  Jp  jz^o  j  c‘jno 

r  +  +  -j2irk,T 

Jdr o  P  (rQ) e  (4.10) 


If  the  density  function,  p (r  ) ,  is  assured  Gaussian  with  a 
standard  deviation  p  ,  then  the  averaged  coefficients  become 


.n 


|R(k) 


1  2  ,,  2  .  2. 
-Zap  (k  +  k  ) 
*o  x  y 


(k)e 


J  n ' 


!  ~  2  ]  ’ 

2 71  p  \  Vkv  d 


(4.11) 


The  Gaussian  assumption  on  rQ  is  reasonable  because  rQ 
often  represents  a  pointing  error  from  some  designated  tar¬ 
get  position.  Then  po  describes  the  accuracy  of  the  pointing 
mechanism.  Substituting  (4.S)  into  (3.42),  the  filter  func¬ 
tion  for  the  threshold  detector  becomes 


q  (t) 


1 


TTf  1 

1  +  ~  K-p 


/T 


(4.12) 
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Rewriting  the  equivalent  signal -to-noise  ratio  in  terms 


of  equations  (4.8)  and  (4.12)  one  lias 


M 

I  a  /r/R 


/hi/  2 


r 

k 


h 


_ 1 

UTK-Bk 
1  ~W2~ 


(4.13) 


It  will  be  shown  in  Part  V  that  may  be  represented  as 


«k  '  °B2n\ 

2 

where  is  the  amplitude  of  the  background  noise  correla- 

2 

tion  function,  is  the  instantaneous  f ield-of- view  of  the 
detector  and  Qk  is  an  integral  function  of  the  optical  and 
detector  parameters  and  correlation  length  of  the  background 
radiance . 

Define  the  following  terms: 

*  S  J.  VT 

Signal  -  to -noise  ratio  (SN’R)  = - — -  (4.14a) 

►/N>  2 


Og0K/T 

Background- to-noise  ratio  (BN’R)  =  —  -  - 

Sn/2 


Using  (4.14a)  and  (4.14b),  d  n.ay  be  rewritten  as 


SNR  AT 


I", 


1  ♦  MQ,  PNR"  I 

*  _J 


(4.14b) 


(4.15) 


Noise  equivalent  irvadiancc  (NF.T)  may  be  defined  at  the 


detector-preamplifier  output  to  be 


(4.16) 


K'NEI  -  J*  •  2Bk,  '■  ATVW 


where  B^  is  the  detector-preamplifier  bandwidth. 


The  SNR  and  BNR  may  now  he  written  in  terms  oi  the  NFI  as 

/R  SNR  «=  “~n *  SNR'  /FTTTT 
Nh  I 

where  NTB  is  number- time -bandwidth  product  and 


Or  _ 

/R  BNR  -  jjyfy-  /RTF  -  BNR'  /RTF. 
Expressing  d  in  terms  of  NEI  one  has 


d  -  SNR'  /RTF  E  | H,  | 
k  k 


1  +  NTB  BNR '  *"  | 


(4.17) 
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V.  SOLUTION  FOR  RFCTANGULAP.  DFTFCTOR 


The  equations  in  Part  II  for  the  detector-preamplifier 
output  voltage  and  the  correlation  function  of  the  background 
must  be  solved  to  obtain  the  filter  function  and  specify  the 
detection  system  performance. 

The  two  general  equations  obtained  in  Part  II  are 


where  H  (r  )  is  the  nth  Fourier  series  coefficient  for  the 
radiation  on  the  spatial  filter  and 


B 


n 


«= 


I  Fo  CV) !  2  WB(k)  J 


2 

n 


d"k 


where  is  the  n*^1  coefficient  of  the  background  correlation 
function  at  the  spatial  filter. 

A  rectangular  shaped  detector  was  chosen  for  the  calcu¬ 
lations  because  it  exhibited  the  simplest  form  of  the 
equations.  For  circular  nutation  the  filter  is  represented 
graphically  in  Figure  5.1. 


Figure  5.1.  Nutation  with  Rectangular  Detector. 


*Ar) 


’  e  20 


T~ 


2no ' 


Its  transform  is 


F0(k)  -  c 


2tt2o2  (k 2  +  k 2 ) 


CS.2) 


The  correlation  of  the  background  radiance  is  assumed  to 
be  of  the  form 


4B(x,y)  -  cj; 


(5.3a) 


where  a  *  and  f?  *  arc  the  correlation  lengths  in  the  x-  and 
y-  directions,  respectively.  The  Kicncr  spectrum  is  the 
Fourier  transform  of 


2  2a  2  B 

“  °b  ~i - ? 


a‘  *  (2mk  )  r  ♦  ( 2  tt  k  )  “ 
*  / 


(5.3b) 


Substituting  (5.1),  (5.2)  and  (5.3b)  into  (2.10),  the 
radiance  from  a  point  target  becomes 


Hn(ro’  *  3  /«  *  X  y  y  e  x  y 


m/e 


j  2lr  (kxxo+kvyo^  sin’Tkxw  sinrkyh 

' e  mk  mV 

x  y 


-1 


-in  tan  (k  / k  ) 


y/kx’  Jn^V^ky)dVV 


(5.4) 


Letting 


ux  .  /IroVx 


l,y  “ 
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equation  (5.4)  becomes 


Vv>  ■  h  °  y 


Qrx’xo^ux*  ^ry  ’yo^uy3 


-jn  tan'  (u  /u  )  // y  /-? y\ 

e  J»[i_  JU*UY  Jduxd"y  t5’5) 

If  rQ  is  random,  the  radiance  from  a  target  may  be 
averaged  with  respect  to  the  density  function  p  (r  )  as  shown 
in  Part  IV.  In  particular,  if  rQ  is  Gaussian  distributed, 
with  standard  deviation  pQ  the  averaged  coefficients  become 


'n  ' 


where 


sirnik  w  sin~k  h  -jn  tan  ( k  / k  1 
- ^  6 


Jn  2TTp7kx  +  kv  dk.xdkv 


(5.6) 


Letting 


2  2  .  2 
Y  ■  pQ  ♦  a 


ux  ■  ^kx 


u  '  *=  /Iryk 

y  y 


equation  (5.6)  becomes 
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H 


n 


'2  '2 


-jn  tan*1 (uj/u') 
e  *  7 


•Jn  jduxduy  C5-7) 

And  finally  substituting  equations  (5.1),  (5.2)  and 
(5.3b)  into  equation  (2.22),  the  correlation  coefficients 
become 


'n 


/• 


4m2c2(k2-k2)  sin2  (irk  w)  sin1,  (irkvjrh) 

e  x  y  - j — 


(^kx) 


(*ky) 


2  2a 

°R  T 


7  7T 


a  +  (27T^'x^  B  +  C^ffky ) 


•> 

-X-  J“ 

2  n 


2  np 


7^ 


Letting 


y|dVky 

r  c  oi 
v  •  «  ) 


Ux  -  2^kx 


U”  *  2rrpky 


equation  (5,8)  becomes 


2  2  2 
oRn  aBo  r  -u  -u 

K  ■  ■  J'  y 


"2  ,sin  it  2/sin  u”\2 

*To 


W 


/To 


u 


X  ) 


/To 


/To 


u 


y  > 


l 


(ao)  2  +  u*x  (Bo  r*uy 


1  t2  p  /"2 

t — rrv  J  -  /u  ♦u 

,2 .  .'.2  n  o  \/  x 


y 


du’’du''  (5.9) 
a  y 


where 


n2  -  w2h2. 


(5.9) 
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''•** *r«v«Pfr4PJUi 


The  integrals  in  (5.5),  (5.7)  and  (5.9)  were  difficult 
to  evaluate.  Although  the  components  of  the  integrand  are 
relatively  simple  functions,  the  argument  of  the  Bessel 
function  prohibited  separating  the  integral  into  two  one¬ 
dimensional  integrals.  No  closed  form  solutions  were  found 
so  the  integrals  were  evaluated  numerically.  The  numerical 
method  used  is  based  on  work  by  Pierce  [9]  who  applied 
Gaussian  quadrature  formulas  to  two  dimensional  integration 
by  integrating  over  a  planar  annulus  in  the  (x,y)  plane. 

Gaussian  quadrature  formulas  are  a  means  of  evaluating 
the  integral  by  summing  weighted  values  of  the  integrand  at 
specific  points.  Fcr  the  one  dimensional  case, 

b  b 

J  g  (x)dx  «  /  w(x)f  (x)ux  (5.10) 

a  a 

where  g(x)  is  the  integrand  to  be  integrated  and  w(x)  is  a 
weight  function  for  which  the  specific  integration  was 
derived.  For  well  behaved  funztions  the  integral  may  be 
evaluated  as 

b  m 

/  g (x) dx  «  A  . f (x . )  ♦  error  (5.11) 

a  i-1  1  1 

where  the  x^'s  are  the  m  zeros  of  the  mt^1  polynomial 

m 

P  (x)  -  n  (x-x.) 

*  i-1  1 

of  the  set  of  polynomials  mutually  orthogonal  on  the  inter¬ 
val  t a  , b ]  with  respect  to  the  weight  function  w(x).  The 
weights  arc  given  by 
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(5.13) 


b 

■  /  w(x)  L^(x)dx 

a 

where  1^  (x)  -  Pm (x) / ( (x • x ^) (x^ ) ]  is  the  Lagrange  inter¬ 
polation  coefficient. 

■Por  a  weight  function  w(x)  and  interval  [a,b]  the  poly¬ 
nomial  Pm(x)  and  its  zeros  xi  and  weight  factors  need  be 
computed  only  once.  For  a  number  of  weight  functions  and 
intervals  the  set  of  orthogonal  polynomials  is  known,  Stroud 
and  Secrest  [13]  give  x^'s  and  A^'s  for  a  variety  of  weight 
functions  and  internals.  The  degree  of  the  formula  deter¬ 
mines  the  number  of  x^'s  and  A^'s.  A  2M-1  degree  formula 
will  have  M  points  and  is  exact  for  polynomial  integrands  of 
degree  2M-1  or  less. 

Pierce  applied  Gaussian  quadrature  integration  to  two 
dimensions  where  the  solution  is  of  the  form 

//g(x,y)dxdv  =  r  I  ^ j £ (xi j  .Yi j ) 

The  integration  is  over  an  annulus  in  the  x,y  plane  with 
inner  radius  R  and  outer  radius  1.  Rewriting  the  integral 
in  polar  coordinates 

1  2  TT 

/  /  rG (r , 9) d9dr  =  l  l  P..G(r,,9.)  (5.14) 

R  o  i  j  1J  J  1 

Pierce  showed  the  summation  could  be  rewritten  as 

4(m+l)  m+1 

l  I  A-B.G(r.  cos6.,r.  sinO.)  .  (5.15) 

i  *=  1  j  =  l  1  J  1  J  1  J 
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where  and  Rj  arc  the  weights  for  the  radial  and  angular 
directions  respectively  ard  4M+3  is  the  degree  of  accuracy. 

For  an  arbitrary  annulus  the  formula  may  be  obtained  by 
rewriting  equation  (5.14) 


2  tt  t2 

I  ■  /  /  f(r,0)  rdrdO 

0  rl 


(5.16) 


and  letting 


p  *  ~ 

r2 


,  2  TT  1 

!  I 

°  rl/r. 


pf(r2p,0)  dpdG 


(5.17) 


The  approximate  solution  can  be  derived  as 
4(m+l)  m+1  j  n - 7-^- 


(5.18) 


where 


■  2ir/4(m+l) 


2  2 

^  ■  wj(r2-r1)/2,  u.  -  weight 


M+l  is  the  order  of  the  orthogonal  polynomial,  in  this  ease 
the  Legendre  polynomial  on  the  interval  (0,1)  andCj  's  are 
the  zeros  of  this  polynomial.  This  type  of  formula  is  known 
as  a  spherical  product  Gauss  formula  [12]. 


I 


of  the  integrands.  One  additional  annulus  inside  the  first 

zero  was  found  necessary  to  account  for  the  — term  when 

a+x 

a  is  small  and  x  (and  y)  approach  zero.  The  distances 
between  annuli  outside  the  first  five  to  eight  zeros  were 
found  not  to  be  critical  and  thus  spaced  arbitrarily. 

The  Gaussian  Quadrature  formula  used  was  a  Gauss-Lcgendre 
24 • point . formula  as  listed  in  [13]  thus  giving  a  degree  of 
accuracy  of  99.  The  integration  is  limited  to  only  the  first 
quadrant  as  mentioned  below.  Thus  for  a  24-point  formula, 
(M+l)  or  576  points  per  annulus  were  used  to  evaluate  the 
integral . 

The  expression  for  6n  is  easily  seen  to  be  an  even  func¬ 
tion  both  radially  and  about  both  planer  axes,  therefore, 
only  integration  over  the  first  quadrant  was  necessary. 

With  a  little  more  difficulty,  f'n(ro)  and  l'n  can  be  shown 
(see  Appendix  F)  to  exhibit  the  sane  property  provided  one 
expression  is  used  for  even  n  and  another  for  odd  n.  3y 
calculating  only  over  the  first  quadrant  computer  time  was 
considerably  shortened. 

In  order  to  solve  for  coefficients  of  arbitrary  order, 
the  Bessel  functions  of  that  order  must  first  be  obtained. 

The  available  library  subroutines  were  found  inadequate  for 
this  use  for  two  reasons.  First,  these  routines  have  an 
upper  limit  of  100  on  the  order  of  the  Bessel  function  th3t 
can  be  computed  and  second,  the  routine  had  to  be  called  for 
each  individual  order. 
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The  recursive  relation 


\(y)  ■  f\tlM  ■ 


(5  19) 


provides  a  basis  for  generating  an  array  of  Bessel  functions 
of  various  orders.  However  for  accuracy  Jn  +  ^(Y)  and 
must  be  known. 

The  most  accurate  method  to  numerically  calculate  Besse1 
functions  of  various  orders  and  arguments  was  found  to  lie  a 
uniform  asymptotic  expansion  involving  Airy  functions  [1]: 


11  ' A£  (n 2 7 ^ A )  -  ak(A) 

— it? —  v£n  ~TfT 
n  kc0  n 


A!(n2/3A)  -  b.(A) 

+  — J7T—  ,.L  — 


n 


k*=0  n 


(5.20) 


where  A^  and  A j  are  two  types  of  Airy  functions  and 


y  a3/2  "  ln 


1  + 


A 


'/TT1 


2  <  1 


aQ(A)  1 


b_  (A)  -  -5/48  A2  ♦  — 
0  SK 


i  rrr 

24  Cl - “ 


- r- 

8(1*4")^ 


z  <  1 


ak(A)  <<  aQ  (A) 

bk(A)  <<  b Q  ( A ) 


k  >  1 


(5.21) 


The  Airy  functions  are  calculated  using  the  equations 


A|  (x) 


~  e"6f(‘i) 


?.  nd 


A  •  (x)  =  -  j  x  1  c'6g(-n 
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where 


6  «  |  x3/2  (5.22) 

and 

x  =  n  A 

The  functions  f(-6)  and  g(-<5)  are  tabulated  and  linear 
interpolation  is  used  to  determine  their  values.  A  is  guar¬ 
anteed  to  always  be  positive  if  n  is  selected  to  be  larger 
than  the  argument. 

By  picking  an  order  much  higher  than  the  argument,  v, 
and  J  using  the  asymptotic  expansion,  the  recursive 
relationship  (5.19)  may  be  used  to  generate  an  array  of 
successive  orders  of  the  Bessel  function  down  to  and  includ¬ 
ing  JQ .  However,  error  builds  up  rapidly  using  this  method. 

A  second  relation 

1  =  JQ(y)  ♦  2J2(y)  *  2J4(y)  ♦  ...  (5.23) 

may  be  used  to  generate  a  normalizing  factor 

C  -  l/[j0(y)  +  2J2(y)  *  2J4(y)  +  .  .  .]  (5.24) 

Multiplying  each  Bessel  function  value  generated  in  the 
recursive  relation  by  C  resulted  in  very  accurate  values 
being  obtained.  The  computer  program  used  to  solve  the 
equations  is  listed  as  Appendix  H. 

To  evaluate  the  numerical  results  a  check  was  used  which 
summed  the  coefficients  anti  compr.cd  this  summation  to  an 
analytic  expression  for  a  sum  that  could  easily  he  solved 
numeri cal lv . 
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For  a  particular  phase  angle  of  the  nutation  cycle,  uQt , 
the  radiance  function  for  a  target  is 


*  £  Hk(r0)c 


jku)0t 


(5.25) 


For  p  point  target  and  a  rectangular,  circularly  nutat¬ 
ing  detector 

-2Tr2u2k2  single  j  2tt  k  (pcosio  t  -  (r  -x  )) 
H(«0t)  «  fe  X  e  X  0  x  0  dk} 


fe 


-2Tt2o2k2  simrhk  j  2irk  (psinw  t  -  (r  -yQ)) 
y  _ _ _ l  p  *  ' 


— C -  e  '  '  dk 

y 


(5.26) 


which  has  the  solution  (see  Appendix  F) 


H(V)  =  i- 


Pl(rx’X0,W)  *  (rv*Xn»'W) 


1  '•‘x  ,Ao« 


-  F2fry,yo’h) 


E2(rv,yoh) 


where 


r  r-  z  -ocosu  t 


E ,  (  r  ,  2  ,  a )  =  Erf 


isi*  -a-1 

/To  2  /To  J 


E9(r,z,a)  =  Frf 


r-z-psinw  t 

_ 2_  +  — — 


/To 


2 /To  J 


and  Erf  is  the  error  function 


2  X  ■  U 2 

F.rf(x)  =  — -  /  e  du  . 
/tT  0 


(5.27) 


By  evaluating  the  analytic  expression  at  the  phase  in' 


stant  (u  t)  that  the  detector  crosses  the  target  and  comparing 
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it  against  equation  (5.27)  for  ar  appreciably  high  n,  a 
reasonable  chock  was  made  on  the  accuracy  of  the  individual 
signal  coefficients. 

For  the  coefficients  averaged  over  r  ,  the  check  was 

made  at  u  t  «  0  where 
o 

CO 

H(o)  -  E  H.  (5.28) 

k  K 

and  the  analytic  expression  is 

H(o)  -  ‘  El(rx'xo>-w)][E3(Vyoh) 

‘  E3Cry -yo-hi] 

where 

E3(r,z,a)  =  F2(r,z,a)l  tc0  (5.29) 

o 

The  rms  background  voltage  was  used  to  check  the  coef¬ 
ficients  in  the  correlation  function.  The  mean  square  volt¬ 
age  is 

00 

eg  «  E{B2(t)}  -  RB(0)  =  E  3^  .  (5.30) 

Appendix  B  shows  that  the  correlation  function  may  be  written 
as 

Rg (t )  -  / | t  (k) | “  lf£(k)  Jo  4-ek  sin  d^k  (16B) 

Substituting  (5,1),  (5.2),  and  (5.36)  in  (16B)  and  evaluat¬ 
ing  at  t  =  0 
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■2  ,  sin*Cnkxwl  c-4’2°2kx 

c» '  0*2a ' — JZ?  ^tt^T 


7  dhx 


.  2  ,  ,  .v  -4/c  k  , 
sin  (rt  h)  e  y 

•o-26  / - - — V--  —5 - j  ^ 

B  (mlO2  BZ  ♦  f 2irk  ) 

/■  y 


(5.31) 


It  can  be  shown  (see  Appendix  F)  that  an  analytic  expres 
sion  is 

eZ  =  4aBOgR  ( w ; a)  R (h  *,  6) 


where 


R(u ,a) 


-  u 

7 


1  -  Erfc 


2v 

a 


2  2 
a  a 


2a‘ 


eay  Erfc 


2  a 
/?y 


ac 


1  -  e 


-^2 


2o 


♦  e*ay  Erfc  |aa  -  ^-j  -  2Frfc(ac) 
and  Erfc  is  the  complimentary  error  function 

Eric  (x)  =  —  /  e'u  du  . 

/if  x 


,1 


(s.:  .) 


(5.33) 


The  computer  program  used  to  solve  equations  (5.6)  , 
(5.7)  and  (5.9)  and  calculate  the  infinite  summations  is 
listed  as  Appendix  H. 
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VI.  NUMERICAL  RESULTS 


Numerical  results  constitute  solving  for  the  signal  and 
background  coefficients  and  using  these  coefficients  to 
specify  the  form  of  the  filter  to  he  used  in  the  threshold 
detector  as  well  as  to  calculate  the  probability  of  detection 
as  a  function  of  s ignal - to -noise  ratio  (SNR),  background  -  to - 
noise  ratio  (BN'R)  and  number  of  nutations  (M)  used  to  make 
the  decision.  The  probability  of  detection  is  determined 
for  the  averaged  equivalent  signal -to-noise  ratio  and  for 
specific  point  targets. 

Many  different  values  for  system  parameters  and  target 

locations  were  used  to  compute  the  coefficients.  A  typical 

set  of  parameters  is1: 

nutation  radius  =  15 
width  of  detector  =  50 
height  of  detector  c  1 
blur  circle  standard  deviation  *  .5 
background  correlation  length  x-direction  =  CO 
background  correlation  length  y-dircction  *  CO 
position  of  detector  in  detector  coordinate  system 
x  direction  =  15 
y  direction  =  0 

target  coordinates  (if  desired) 
x  direction  c  0 
y  direction  »  0 

pointing  error  standard  deviation  =  7.5 
A  detector  utilizing  the  parameters  above  would  trace  out 
the  area  in  the  image  plane  shown  in  Figure  6.1. 


1 


All 


units  arc 


in 


r.  i  1  1  i  r  a  d  i  a  n  s  . 
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TABLE  6.1 


Calculated  Sum 

Actual  Sum 

Signal 

(point  target) 

.6826894 
.6826895 
.  5438002 

.6826993 
.6826994 
.  5438087 

Signal 

(averaged) 

.2990664D-6 
.2328138D-4 
.  2651  595D- 1 

.  2991232D-6 
.2328142D-4 
.  26  51  S9SP-1 

Background 

.  601  34  77 
.2347420 
.8335525 

.6012230 

.2347420 

.8319764 

Besides  the  general  shapes  of  the  background  and  signal 
spectra,  variations  in  the  spectra  due  to  detector  size, 
background  correlation  lengths  and  coordinates  of  the  detec¬ 
tor  are  of  interest.  Curves  in  the  following  figures  show 
the  effects  of  these  variations. 

Figures  6.2  through  6.5  show  the  envelope  of  the  one 
side  of  the  two-sided  background  power  spectral  density 


sbM 


oo 


*=  l  6,  6  (w-kuj  ) 

k  k  0 


where  is  normalized  to  2r .  The  curves  of  Figures  6.2  and 
6.3  illustrate  the  effect  of  varying  the  background  correla¬ 
tion  lengths.  In  Figure  6.2  the  correlation  lengths  are  the 
same  in  both  x-  and  y-directions .  Practical  measurements, 
however,  indicate  that  correlation  lengths  in  the  horizontal 
(x-)  direction  may  be  longer  than  in  the  vertical  (>’-)  direc 
tion.  The  curves  in  Figure  6.3  arc  for  this  case. 


5  8 


The  curves  of  Figure  6.4  and  6.S  relate  the  spectrum  to 


changes  in  size  of  the  detector.  In  Figure  6.4,  only  the 
detector  with  RNR  equal  to  9  and  number  of  nutations  equal 
to  100  are  shown  in  Figures  6.11  through  6.13.  Each  figure 
shows  the  envelope  of  one  side  of  a  double  sided  spectrum. 

Two  cases  are  considered.  First,  threshold  detection  systems 
designed  for  different  Gaussian  pointing  errors  but  with  the 
same  size  detector  are  shown  in  Figures  6.11  and  6.12  and 
second,  detection  systems  designed  for  the  same  printing 
error  but  utilizing  different  size  detectors  arc  shown  in 
Figures  6.12  and  6.13. 

The  curves  shown  in  Figures  6.14  through  6.16  show  the 
matched  filter  to  a  particular  point  target  located  within 
the  scan.  Notice  each  of  the  filters  exhibits  a  band  pass 
characteristic.  This  is  because  of  the  additional  noise 
component  (background).  The  oscillations  observed  in  the 
point  target  spectra  of  Figures  6.6  to  6.S  have  been  sup¬ 
pressed  to  provide  a  look  at  realizable  filters. 

The  probability  of  detection  for  the  threshold  detector 
was  shown  in  Part  IV  to  be 

l  |f  |2  - T 

k  *  1  ♦  MQkBNir 

Figures  6.17  through  6.20  show  as  a  function  of  SNR, 

BNR  and  M. 

The  probabilities  of  detection  for  point  targets  located 
at  specific  points  in  the  scan  have  also  been  plotted.  The 
probability  of  detection  for  a  specific  point  target  is  shown  as. the 


Qd  -  SNR  /R 


width  of  the  detector  is  varied.  As  width  is  increased,  so 
is  the  area  scanned  by  the  detector.  Figure  6.5  shows  the 
difference  in  spectra  between  a  rectangular  detector  and  a 
square  detector  which  have  approximately  the  same  surface 
area.  It  should  be  obvious  however  that  the  rectangular 
detector  scans  much  more  area  per  nutation  than  a  square 
detector  of  identical  surface  area. 

The  envelopes  of  signal  coefficients  (magnitude  only) 
for  a  point  target  are  plotted  in  Figures  6.6  through  6.8. 
Phase  information  depends  only  on  the  location  of  the  target 
relative  to  the  initial  point  for  the  nutation  cycle  which 
is  along  the  x  axis.  This  has  not  been  plotted.  Changes  in 
the  magnitude  due  to  changes  in  dctectoi  size  are  shown. 

Figures  6.y  and  6.10  show  the  envelopes  of  coefficients 
for  the  signal  that  is  averaged  with  respect  to  a  Gaussian 
pointing  error.  The  bandwidth  of  the  averaged  signal  is 
much  less  than  a  point  target  as  expected.  The  figures  show 
the  spectra  for  two  different  standard  deviations  of  point- 
ing  error. 

Matched  filters  (magnitude  only)  for  the  threshold 


=  SNR'/CT.- 


£  fk^ro^k  i  +  MQj-BNR2 


r  |  f  I  ^  1 

1  I  rk  !  - 7 

k  K  l+MQkBNFT 


where 


SNR 


1 

M 


M 

T, 

n=l 


\ 

2 
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Figures  6/21  through  6/24  show  Qj  as  a  function  of  SNR', 

BNR  and  M. 

The  parameters  of  the  nutating  optical  system  have  been 
varied  to  show  their  effect  on  and  Q^. 


Detector  parameters:  RH0=1C,  SIG=.5,  w=10. 


Common  detector  parameter 
Bar,  '”-ni,nfl  na r ame  t e r r> :  a 


Common  detector  parame 
Background  parameters: 


Dev.  of  Gaussian  Pointing  error 


shold  Det 


Figure  6.17.  Probability  of  Detection  vs.  Signal - to-Noise  Ratio 
Detei.^;  ji  j  rameters  :  RHO-15,  SIG=.5,  w=10,  h  =  l, 

xo=15>  >VC 

Std.  Dev.  of  Gaussian  pointing  error  =  5.0 


obability  of  Detection  vs.  Si gnal - to-Noi se 
f.-rtnr  Parameters:  RIIO'15,  SIG-.S,  w=30, 
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CONCLUSIONS 


The  problem  investigated  in  this  thesis  was  to  determine 
an  optimum  processor  to  be  used  in  detecting  targets  with  an 
infrared  nutating  system.  The  nature  of  the  system  divided 
the  problem  into  two  areas:  optimising  a  spatial  processor 
according  to  size  and  shape  and  optimizing  a  temporal  pro¬ 
cessor  that  takes  the  output  time  varying  voltage  and  decides 
if  a  target  is  actually  present. 

General  equations  for  the  output  from  a  nutating  detector 
were  known  from  Samuclsson's  work,  however  these  equations 
had  not  been  solved  for  a  specific  detector.  At  the  same- 
time  it  was  observed  that  the  fuiiVi  of  the  covariance  function 
for  the  noise  offered  an  easy  solution  to  integral  equation 
of  the  Karhunen- Locvc  expansion  which  made  statistical 
detection  theory  appealing.  Some  work  had  been  done  on  the 
form  of  temporal  processor  using  statistical  detection  theory 
but  this  was  limited  to  only  one  nutation  because  of  the 
common  background  noise  component  between  nutations. 

Harger's  derivations  provided  a  means  for  describing  a 
temporal  processor  using  statistical  detection  theory  that 
based  its  decisions  on  multiple  observations  and  included 
the  background  noise.  His  work  was  extended  here  to  include 
unknown  parameters,  amplitude  and  position,  which  more 
closely  characterized  the  system.  This  extension,  called  a 
threshold  detector,  was  shown  to  be  optimum  in  the  case  of 


>y*>KTiin  'ii'wiwl.u 


small  s  1  gna  1  - 1  o  -  no  i  s  e  ratios  where*  detection  is  most  diffi¬ 
cult.  The  derivation  also  used  the  Gaussian  assumption  to 
describe  the  background.  This  may  or  may  not  be  correct. 

But  in  practice,  the  detector's  performance  may  not  suffer 
greatly  if  this  assumption  is  wrong. 

To  implement  the  threshold  detector,  the  integrals 
describing  the  dctecor  output  were  solved  using  the  Gaussian 
quadrature  method  of  numerical  integration.  Checking  the 
computed  coefficients  by  a  summation  provided  a  means  of 
accepting  the  validity  of  the  integrations.  The  form  of  the 
threshold  detection  system  using  a  rectangular  detector  was 
determined  and  the  frequency  spectrum  of  the  opt  in urn  filter 
was  shown.  To  specify  performance,  the  probability  of 
detection  is  plotted  against  signal -to-noisc  and  background- 
to-noise  ratios  and  the  number  of  nutations  on  which  the 
decision  was  based.  The  signal  spectrin  for  a  point  target 
and  background  power  spectral  density  is  also  plotted. 

A  designer  may  find  these  result.',  useful  ir.  developing  a 
system  or  investigating  the  performance  of  an  actual  system  s 
compared  to  the  optimum.  He  may  also  use  these  results  to 
determine  the  sire  of  a  rectangular  spatial  filter  to  be  used 

Future  work  should  include  extending  these  results  to 
spatial  filter.--  that  are  circular  or  elliptic.  The  same 
numerical  integration  algorithm  could  probably  be  used.  The 
threshold  detector  was  found  to  be  highly  sensitive  to  back¬ 
ground  correlation  lengths.  Thus,  some  form  of  adaptive 
processor  should  be  included  to  minimise  the  effect  cf 
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background.  At  the  present  time,  a  measurement  program  is 
underway  at  Naval  Weapons  Center,  China  Lake,  California 
to  determine  average  background  correlation  lengths  but  for 
different  environments,  these  correlation  lengths  could  be 
expected  to  vary  considerably. 
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APPENDIX  A 


DERIVATION  OF  SIGNAL  SHTTRAI.  GOFF  1: 1  C I  I'N'TS 


The  radiance  on  the  detector  is  periodic  because  of  the 
nutation  therefore  it  can  he  expressed  in  a  Fourier  series, 

00  jnu  t 

H(t)  «=  I  H  e  0  (1A) 

n 

vherc 


Vt  /  ,l(t)c 


But  , 

H(t)  =  /n' (r)T(r-p(t))  d2r  .  (3A) 

Substituting  (3 A)  in  (2 A) 

f  ->  if*1'  -Jnu.0t 

Hn  =  Jd"r  N '  ( r )  ±  J  dt  'Cr-c(t)  e  (4A) 

The  transmittance,  t,  nay  be  expressed  in  terns  of  the 
inverse  Fourier  transform,  assuming  an  infinite  image  plane, 

r  ,  -j  2-rrk •  fr-p  (t)) 

TCr-o(t))  -  J  d-k  r*fk)e  -  ~  .  (5A) 

Rewriting  the  exponent, 

k*(r-p(t))  =  k-r  -  pk'xcos.ct  -  ol'vsinuot  (bA) 

Equation  i-h\)  can  be  rewritten, 
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where  background  covariance  function  on  the 

image  plane. 

It  is  now  convenient  to  express  the  integrand  in  terms 
of  two-dimensional  spatial  Fourier  transform.  To  begin  with 
,  j  2rk • r  9 

t (r )  -  j T(k)e  ~  -  dzk  (6B) 

where  t (k)  is  the  Fourier  transform  of  the  aperture  function 

t(rj  and  k  ~  (k  ,k  ).  Then 

-  -  x  y 

T  [rre(^+T)]  t  [r2'£ 

?  f  ,  j2n(k, •r1-k7-r7) 

=  Jd  kxT (kj)  Jd-k2T*(k2)e  -1  ~2  ~2 

-  j  2  r  fk .  0  (tri  ) -k2C  (-t)l 
•  e  L-i  *  J  ( 7  B) 

Substituting  (SB)  and  (7B)  in  (4B)  and  interchanging  the 
order  of  integration, 

Rb(t)  =  /d2kx  tfkj)  / dZk2  T*(k,) 

i  rT  -J2r[k  -P(t*T) -k9-o(t)l 

.  i  dt  e  L  1  '  -  J 


•  fd2r. 


j  2rbi  ‘£i 


/d2r 2 


■  j  2irk0  -  r 


The  last  integral  becomes,  with  R  =  r.-r. 

-  -  i  —  i 


.  ,  ‘  j  2  r.  k , '  R~|  -j2-k9-r 

d“Rt  '  (R)  e  ~  e 


But  the  bracket  above  is  the  Fourier  transform  of  g(r),  or 
the  Wiener  spectrum  of  background  on  the  image  plane,  Let 


-  j  2  r  k  ,  •  R 

W’(k:)  -  Jl’( R)e  --  ~  d-R 


(10B) 


9 


If  FQ(k)  is  the  Fourier  transrorm  of  the  point  spread  func¬ 
tion  and  Ku(k)  is  t  lie  V.'iencr  spectrum  of  background  on  the 

D  - 

object  plane,  then 

K'(k)  =  |F0(k)|2  K. (k)  (11B) 

This  relation  is  analogous  to  the  output  noise  spectr  i 
expressior.  du  to  a  noise  input  to  a  linear  system.  In  terms 
of  Wjj(k)  the  last  line  of  (SB)  is 

r  j  2  n  (k,  -k7)  T,  ? 

h’g  (k)  J c  -1  ^  d ir1  -  Wy  C^)  <5  (kp-^)  (HB) 

Substituting  this  into  (8B)  and  performing  the  integration 
with  respect  to  k?, 

rb(t)  =  /d2k  | t  ( k )  |  2  W'(k) 

i  f  T  j  Zrrk-  [C  (t+O  -o  (t)  ] 

•  y  J  e  ~  dt  (13B) 

o 

Since  p(t)  =  (p  cos^t,  psinc0t), 

k  •  [c  (t+t) -p  (.t)  ]  =  kxccoswQ  (t  +  -  j  +  k  cs  inuQ(t  +  7) 

-kxpcos^Qt  -  ky.sina'Qt 

-  A(T)cos^Qt  +  B  (  t  )  s  inajQt  =  c  (  t  )  c  os  (oQ  t  -  „  ( t )) 

whei  e 

Aft)  =  p  [  k  cost-  t  +  k  sin-  t  -  k  ] 

^  x  o  y  n  x 

B(t)  =  c[-kxsinuoT  +  k  c o s c o t  -  k  ] 

c  ( o  *  v^(0  +  k:(T)  =  P  N,f2(l  -  cos,,oO  (X,  -  k2) 


Therefore  the  last  integral  of  ( 1 3 B )  is 


,  -T  j2^c  cosf^T  +  'J') 
y  J  e  dt  c  Jq(2ttc) 

1  o 


(15B) 


where  J  '  • )  is  the  zeroth  order  Bessel  function.  Therefore, 

V.  ' 

the  correlation  function  Rg(t)  is  given  by 


Rb(t)  =  /|t(k)!2  K’(k)  J 0 f 4 n pj<  sin  ^d2k  (16B) 


where 


k  =  J k2  +  k2  and  d2k  =  dk  dk  . 


Finn  .ly,  the  power  spectrum  coefficient  S  is,  from  (IB) 


en  =  / I t (k) | 2  (k)d2k 


w  r  \  •  j  nu;  t 

2a  sin  —  I  e  0  dx 


(17B) 


where  a  =  2^kp.  The  last  integral  becomes  with  x  =  u)qt/2  , 


1  f" 

—  j  JQ(2a  sin  x)  cosZnx  dx 


1  f* 

j  -  I  j  (2a  sin  x)sin  2nir  dx 

J  IT  J  O 


(18B) 


The  cosine  integral  equals  J^(a)  and  the  sine  integral 


vanishes  [2].  Therefore  on  using  (18B) 


/i  T  (k)  |  2  |F0(k)|2  ''  B  ( ~  )  ^  n  (  2 11 P— )  d  2  ^  • 


(19B) 


e 


n 


This  is  the  general  power  spectrum  expression  for  a  nutating 

detector  output.  The  power  at  frequency  is  expressed  in 

terms  of  aperture  transform  r(k),  optical  transfer  function 

Fo(k),  Wiener  spectrum  of  background  W^(k)  and  nutation 
0 

rad i us  through  J^(2upk).  The  integration  is  over  the  entile 
k-pl ane . 


APH1NDIX  c 


KXPtCTATION’  OP  A12 


The  essential  calculation  required  to  find 

A12({Zn^-’  =  rB[A12(Un}'R)]  (1C) 


for  Gaussian  background  B  is  the  expectation 


■ 

T  T 

/  A(t)B(t)dt  -  c  /  B2(t)dt 

exp 

o  0 

. 

where 


c  =  M/N 


and 


A(t)  -  £  E  (Z  ft)  -  S  ft)) 
N  n=l  n  n 


C2C) 


B(t)  m3y  be  expanded  in  its  Karhunen - Loeve  (K-L)  repre¬ 
sentation  provided  the  mean  square  value 

E(B2(t))  <  -  (3C) 

and  thus  the  integral  equation 


T 

X.<t>-(t)  *  /  RH(t,u)  $.(u)du  0  <  t  <  T  (4C) 

11  0  »  i  -  -  . 


can  be  solved.  are  called  the  eigenvalues  of  the  equation 

and  4^(t)  arc  called  the  eigenfunctions. 

Using  Mercer's  theorem,  the  correlation  function  of  the 
background  nay  be  expressed  as 
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-  Z  4*  *  Cu ) 

k 


(SC) 


The  K - 1.  expansion  of  H  ( t)  is 


B (t)  -  Z  bk(^Kf.t) 

k 


with 


bk  c  /  B(t)  $,*  (t)dt 


0  <  t  <  T 


(60) 


vhcre  the  coefficients,  bj,,  are  uncorrclated  and  because 
B(t)  is  Gaussian  the  coefficients  are  statistically  iiulc 
pendent  with  means  uk  and  variances  X^. 

Also  expanding  Aft-), 


A(t)  =  Z  a^^j.  (t) 

1c 


with 


*  /  A (t)  *£(t)dt 

K  o  K 


0  <  t  <  T 


(7C) 


Rewriting  in  terms  of  the  expansions 


I  =  E, 


exp 


i  d  ll  anb5  *  r \l  W  dt 

o  n  k  n  k 


J  OD  OD 


C  /EE  b  bfc  <  x(t)c-*(t)  dt 
o  n  k 


(8C) 


T 

<5XV  =  /  <fx(t)'i*(t)dt 
o  ’ 


Using 


bn 


exp  jf  anb*  *  i  W„)  ■  C  !•  I.nbj 


Uxinjj  the  i mlcpcmlence  of  the  coefficients 


j  -  n  tibn  exp 


L  a  b*  ♦  7  a*l:  -  c 

2  n  n  l.  nil 


exp 


n 


2  t 


n 1 


4c 


n 


l:b„  cxf 


n 

2 

n 

-c 


V2 


n 


7c 


For  Gaussian  random  variables 

cxp|w  mxm*/(l-2KA) 


vhc  re 


!•;  (wxx*) 


mx  *  F.(x) 


T  r  2wT 


A  -  EKx-r.i  )  (x*-m*)| 

A  *' 

and  w  is  a  constant, 
let 


*  '  bn  '  an/2c 


and 

then 


mx  "  pn  ‘  an/2c 


(9C) 


(IOC) 


(11C) 


Ebn  exp 


-c 


an‘  2 


_  CXP  l'c|  Van/?c^(W5jj)  (12 


(12c) 


(1  +  2cXn) 


Substituting  (12C)  into  (IOC),  and  separating  terms 
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B 

n  (i  ♦  2cx  ) 

n 


*1 


exp 


i  w  !  a  |  2  A 
1  v  '  n '  >i 

y  i  irr- 

n  n 


exp 


*»  u * C ii  -a  / c ) 
c  n '  n  n  1 


1“T-2FT 


exp 


c  :  %(l,n'an/c) 

I  *  “T~TcT— 

n  n 


This  result  can  he  written  in  closed  foist  by  using 
definition  A  *  (a,<J>n) 

,  ~  |al2x  i,  T  T 

1  T'  T"+  "2c~A~  *  *4*M  ^  ^  dt2A^t2^>q^tl  *l: 

n  n  o  o 


whe  re 


qttl,t2^)  n~  1  Tel"  *nl  1 


X 


is  such  that  operating  on  with 

/  /  RB(  W  *  WV 


(•) Jtjdt, 


and  using  Mercer's  thcoren  yicls  • 


♦  ;  VW  q(t3^1)dt5  -  Rgdj 
o 


Likewise , 


E 

n 


wn<>n'an/c)  +  lin  K' an/c) 


1  ♦  2c>. 


n 


1 


/ 


’  U  U<0 


-  ^-Wt)  dt 


(13C) 

tlic 

?)  (14C) 

USC) 

USC) 

,t:)  U"C) 

( 1  8  C ) 
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vihc  rc 


hCO 


00 

l 

n 


2q 


n 


T~+""2cX”~  *n(t) 
n 


(19C) 


such  that 


T 

Nh(t)  ♦  2M  I  Rg(t3>t)  h(t3)dt3  -  Mmb(t)  (20C) 

Taking  the  logarithm  of  both  sides  of  equation  (13C)  ar.d 
using  the  equations  above 


in 


I  (1  +  2cX  ) 

n 


+  m  f  dli 

o 


M 


£  ^1  (W  * 


T 

I  dt. 


I  ‘W  -  V2)) 


q  ( 1  j .  t , ) 


T  1 

mb(t)  *  }\ 

[  ,  M  ]1 

/  1 

0 

[  •'  Jl '  1  J 

h(t)dt 


in  I  «  E  (1  ♦  2bXn) 
n 


-1 


1 


T  1 


*  ®  0  0 


M 


M 


M 


M 


I  Z£(t  )  Z  Z  (t 2)  -  I  s  (t  )  I  2  (t,) 
W  1  m=l  m  4  i«l  *  x  m= 1  m  “ 


M  M  M  M 

=  W  £?.W*  .?1s‘(ti)Bl!:1s»(t2) 

i«b1  Z  =  1  Z  =  i  m=l 


q(tj ,t2)dt1dt2 


T 

*  /  m,  (t)h(t)dt 
o 

T  M 

.  2  Z1  r  (Z. (t)  -  S„ (t))  h(t)dt  .  (21C) 

M  o  i  =  l 
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The  likelihood  ratio  desired  is 


Al?({Zn})  “  cxp 


2 


M  T 


2n(t)Sn(t)dt 


1  M  T  2 
i  T.  f  S‘(i)dt 

N  n-1  o  n 


*  I 


(22C) 


Taking  the  logarithm  of  both  sides  and  substituting  the 
expression  for  It  (I), 


M  T 


A12«V>  "  tnlx  ' 


1  M  T  2 
"  /  S‘(t)dt 

N  n=l  o  n 


+  L  tn  (1  ♦  ~  \.)'1 
k  N  K 


,  T  T  ,  M  M 


>1 


M 


M 


M 


1;  W  Z  W  '  1  W  1  Z* 1  ^ 

it  ”  1  1  m«l  *  m-1  m  c  t-1  x  A 


M  M  . 

*  Z  S  o  ( t  -i )  E  S  (td  <l(t.  ,t7)dt,dt, 
£-1  1  m-1  m  1  1  1  -  1  k 


T  7  T  M 

-/  m.(t)h(t)dl  ♦  n  /  1  (Z£(t)  -  S£(t))h(t)dt 

o  b  ‘  o  Jt*l 

(23CJ 
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1 

1 

J 


KAKIIUNIiN- l.OIlVti  UmUiSl.NTATI  0\  01  1-QI1I VAIJ-NT 
- S 1  GN'A’L  -  T0~-~WsF7aTf7) — - 


The  variance  of  tlie  threshold  statistic 


M  fT 

W  -  l  ]  Z  (t)q(t)dt 
n*l  o  11 


where 


q(t) 


2 

N 


f(t) 


•/ 


T  . 


f  (u)q(u,t)du 


(ID) 


under  hypotheses  llo  and  1.1  ^  is 

rl 


VAR 


,  mj  q2(t)dt  ♦  M2  / dfq(t)  /duq(u)K  (t,u)  (2 


D) 


Substituting  (ID)  into  the  first  term  of  the  variance 
expression 


«)it 


2M 

N 


r  ^ 


T 


f2(t)  -  2  f  ( t)  /  f (u)q  (t ,  i  )du 


♦/  /  f (U) f (v)q (t,u)q(t ,v)dudv 

O  0 


(3D) 


Using  the  Karhuncn-Loeve  representation 


4 

-3 


,  ,  q"(t)dt  *  I  I  fkf  *  /  «k(t).t*(t)dt 

2  -  *  k  t  K  *  O  K  *■ 

•  «  -  .  A  rT  -T 

-2  l  ]  hn  7--m Vj  i 

k  x.  1  +  -v-r-  A  o  o 

Is  n 


iu 
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Co  ci'  c*-'  ©v 


2m 


2  m 


*  r.  r.  i;  ?.  f.f* 


4.  in  ^  L.  Ill  •> 

*fT  Am  TT  An 


k  i  m  n  k  *  1  *f  \x  1  ♦  f^n  o 


K—  / 


Using 


/  *k(»HrVl0<H  /  <^,4  (v)  <•*  (v)dt 


(41)) 


‘nk  '  /  ♦11(t)*J(t)Jt. 


then 


t  /  2M  ,  /  2M  . 

x/  q2(Od«  ■  $  r.  fkfj-  f  J  ■  2  -X*  *  '  T  k 


2  ./  *  -  FT  ;  lkxk  \  J  ‘  -  r~m\  . .  + 

°  K  '  1  V  T  Ak, 


2m  r 


“  -FT  J-  VI 


i 

1  +  TT~*k, 


(SD) 


The  second  term  may  be  expressed  in  a  similar  manner  to 


give 


VAR 


1 


00 

2  M  y  r  r  * 

7  ii  177?  » 


*  (x)  H  fkfk(— :V~  1  u 

"k/  '  k  V1  X  'k. 


VAR  ■  X  t  Vk 


1 


k  1  4  xr  A 


TM 

TT  Ak 


(6U) 


The  equivalent  signal  to  noise  ratio  in  Part  II  is 


M 

E  a  '  f(t)q(t)dt 
n  =  1  n  o 


r1 


^  VAR 

The  K-L  representation  of  the  numerator  is 
M  ,-T  -  M 


(  7D) 


l  a  I  f(t)q(t)dt  =  l  a 
n=  1  o  n  =  l 11 


s  £  lfki^(1  +  rr  H 


•  (8U) 
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NON'-Zr.KO  H  UMS  I  OR  NUMl'  V.iCA1.  1  NTlC-llA'i  I  ON 

liquation  (5,S)  for  the  radiance  function  for  a  point 
target  may  be  written  as 

,2  r  1  ?.  -  j  ( i*  Q  ) 

1  n  JL_  f0‘*  *>’  c  X  sy  smax  s i n By 
ti  <  v 


where 


jn<“ 


'Ac 


y-  dxdy 


a  K  W'/ZTo 


and 


p  "  n  /  c  l  q 
c  *  /Ip/o 

P  c  V  (Vxo} 

Q  ‘  T'  (ry'yo} 


-1 


4  *  tan  y/x  (IK) 

Let  E(x,y)  be  that  part  of  the  integral  tnat  is  even  about 
both  planar  axis  which  is 


•  n  2  2 

B Oc.y)  ■  i,  c'x  l‘"a 

TI* 


x  sinSv  T  /  r7  7  \ 

-  -7—  J„  (c  V.X  *y  j.  (21) 


Using  Lulcr's  relation  the  integral  becomes 
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I  =  Jo(x,y) [cos (px+Qy)  +  j  sin(rx«Qy)  1  lcosn'+i  -  jsin  n<t)dxdy  (30) 
Using  trigonometric  identities 


U(x,y)  E  T.  dxdy 

i  =  1  1 

where 


Tj  =  cos  P x  cos  Q y  cos  n  <♦■ 

T ^  =  -sinPx  sin  Qy  cos  nsf 
Tj  =  -sin  Px  cos  Qy  sin  n$ 

=  -cos  Px  sin  Qy  sin  n$ 

Tg  =  -j  sin  Px  cos  Qy  cos  n$ 

=  -i  cos  Px  sin  Qy  cos  nf 

T-,  =  -j  cos  Px  cos  Qy  sin  n<f 

TB  =  *j  sin  Px  sin  Qy  sin  nc  (4E) 

To  save  computation  time  it  is  necessary  to  determine 
which  of  the  vight  integrals  are  always  ;ero  and  for  which  \ 
they  ate  zero. 

Observe  that 


cos  n$  =  R  e^n<^  =  R  (e^) 
e  e 


=  Re(cos4  +  j  sin  4) 


n 


=  R 


\ 

l _ 

2.  2 

x  +  v 


r  2  2 

:x  +  v 


1  - 


y 


■'  *>  0 


\ 


;x  +y  / 


n/ 


Rc  (x  ♦  iy)  . 


(50) 


10  0 


Likewise 


1 

s :  n  n  4>  =  I  — j — 2' 

\xz  +  y 


n/2 


I  (x+jy) 
nr  J  ’ 1 


(6L) 


The  first  term  in  brackets  is  always  even  and  may  be  made 
part  of  E(x,y).  The  second  term  in  brackets  may  be  expanded 
by  the  binomial  expansion  with  the  general  term 


Cx)n-k  (jy)k 


The  first  integral  of  (4k)  becomes 


(7k) 


lj  «  Jk  (x.y)cos  Px  cos  Qy  Re(x  +  jy)n  dxdy 


(8k) 


For  Ij  to  be  non  zero,  Re(x+jv)n  must  be  even  in  x  and  even 
in  y.  Pron  (7k)  this  is  true  only  of  even  n;  for  odd  n,  Ij 
will  always  be  zero. 

Examining  the  other  seven  integrals  yields 

Jb(x, y) [cos  Px  cos  Qy  cos  n* 

+  j  sin  Px  sin  Qy  sin  no]dxdv  n  even 


T  =  < 


|k(x,y)[cos  Px  sin  Qy  sin  no 

+  j  sin  Px  cos  Qy  cos  nfjdxdy  n  odd 

(9E) 

Observing  that  symmetry  about  the  x-  and  v-axis  is  necessary 
for  this  integral  to  be  non-zero  and  integrating  only  when 
the  integral  is  non-zero  permits  one  to  integrate  only  over 
the  first  quadrant. 


107 


APPFNDl X  F 


SUMMATION'  CHLXK 


The  signal  for  a  point  target  was  shown  in  Part  11  to  be 


H(w0t)  c  /n'  (t)t  (r-p(O)  d2r  . 


Assuming  an  infinite  image  plane,  H(u3Qt)  may  be  written  in 
terms  of  the  Fourier  transform  as 


f  -  j  2nk  •  p  (t) 

j  t)  ■  J  N*  (k)T*(k)c  ~  "  cTk 


For  a  rectangular  detector  (2F)  may  be  written  as 


HCucO  =  /  c 


2 11  0  ‘'Kx"Ny 1  sinrwkx  sin~hkv 


{  j  2nk- (p(t) -r  )  -  j2a  (kxx0+ky>'o^  d"k  (3F- 

Jc  •  e  7 


x  c  /?  oTikx 


y  =  /7  cuky 


and  noting  the  integral  may  be  separated  into  twoparts 


2 

,  .  S  1  (  -x  wx 

H(w  t)  =  —  ie  sin  -  e 

o  11  J  r*;- 


j~  x  (pcosm  t  -  (r  -x  ) ) 


o  v  X  O' 


Sin  ST 

,  ,  ,2  /2c  j—  y(psinu  t  -  (r  ->■  )) 


Since  the  integration  is  -*•  to  +°°  only  even  integrands  will 
be  non  zero.  Rewriting  (4 F) 


H(«0t) 


•if 


s  in 


vx 


cosfjp  xlpcos^t  ■  (rx-x0))]dx 


If 

IT  _ 


*  -y2 

e  1 


s  in 


hv 

cosf^|  y  (psinu^t  -  (r„*yo)]j  dy. 


The  general  form  of  the  integral  is 


(5F) 


..2 


R  -  R(B)  =  £/  e'u  •S-1-n^  cos(Cu)  du. 


(6F) 


Differentiating  with  respect  to  B 


— =  ~  j  e’u  cos  By  cos  Cu 


du. 


(7F) 


But  the  solution  to  (71)  is  known  (2]: 


dR(B) 

dB  2/m 


-  (B-C) 


2 


♦  e 


.  (B*C) 
4 


2 


(8f; 


Since  the  constant  of  integration  R C 0 T  -  0, 


R(B) 


2  /r 


/' 


■Cy-c) 


f 


(w  +  c) 


dv  +  J  e 
o 


dw 


(9Fj 


Let 


P  = 


v-c 


and 


1  09 


w  +  c 


then 


R(B) 


1  _L_  f  2 

2  >/v~  J-C/2 


•  P 

c  dP  + 


-hf 

/tP  J-( 


e‘Q  dQ 


1  r  _  ,B-C,  .  r  f  ,B+C\ 

2  1;rf  (-y)  +  Erf  (-=-) 


(10F) 


where 


Erf(z)  =  —  /Z  c'yZ  dy  . 
/rT-  o 


Substituting  the  result  (10F)  into  (SF) 


(1  IF) 


„CV)  =  ^  ( Erf  lff[ 
-Erf  [ec< 


OCOSUot  -  Vxoj  +  v 


pcosoy  -  r  +x 


X01  -  4:  w) 

0j  v‘2a  I 


*Ert  (-"I  '  ry*yol  *  7^:  h) 

-Erf  [rsi„V  -  Vyo]  -  -i-jj.  (12F) 


The  mean  square  voltage  of  the  background  may  also  be 
solved  in  a  like  manner. 

From  Part  V,  the  mean  square  voltage  was  shown  to  be 


*  j ! T (h) ! “  |f0(k)!2  KB(k)  d2k 


which  for  a  rectangular  detector  becomes 


(131) 


-2  , 

=  V2* 


r sin-(rkxK)  c  "  "x 

J  7i  “k1'  a*"  +  (2irk 


- 4r  2o  2k2 


(2nkx) 
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(14F) 


9  ,22.  2 
.•  2,,  -4toK 

sin  (i'k  h)  y 

1 0 u  ‘  28  / - 1 — r-“ - -j - 9  dk 

y 


/s  i  n  2  (  t  I 


Letting 


2  nk 


the  general  form  of  the  integral  to  be  solved  is 

mm  ■  l  f 

o  t 2  /  4  a  ♦  x 

Taking  the  first  and  second  derivatives  one  has 


dR(n )  1_  f°“  s  i  n  ( 2  y  t )  c 


2.2 

-o  t 


du 


IT  J 


a2  +  t2 


dt 


and 


d2Riu)  _  2  r  cos(2ut)  _-o2t2 
'"Jo  02-t2  C 


The  solution  to  the  second  derivative  is  known  [21 


d2R(u)  e°~a" 
du2  ‘  2* 


Erfc  (ao  -  -y-) 

*»  ■J 


*  e 


ay 


[c-“ 

Erfc  (oo  *  jj)] 


where 


?  f°°  -y2 

Erfc (;)  =  —  J  e  y  dy 


/tt  z 


is  the  complimentary  error  function. 

Kith  lengthy  and  tedious  calculations  it  can  be 
by  using  the  relation 


(15F) 


(16D 


( 1  7  F ) 


(18F) 


(191) 


shown 
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APPENDIX  li 


f. 

C  COMPUTER  PROGRAM  TO  CMCULATF  THE  fOU'lER  COE FF  I  C  1  ENTS 

C  no  the  iiackgrodn  j  cokrflt  Hun  funct  it  n.  a  point  target 

f.  ANO  AS  AVFF  A  >..D  TARGET  GVt  w  A  GAUSSi.O'  POINTING  EkkLR 

C  FOW  A  t4lll  AT  1  .\G  OPT  UAL  SYSTiM  W  1 T  M  A  RECTANGULAR 

f.  DETECTOR. 

C 

r. 

c 

r.  AA  =  R ACKGROUNG  CORRELATION  FUNCTION  COEFF (  ZERO  ORDER) 

C  Al  PA  =  i)  AC.KOR  (  i'JNO  CORRELATION  GO  l  l;  F  I C  I  E  N  T  S  (VECTOR) 

C  F  NO  =  Sl(>NAl  lOJ^I  EF.  CGI  FE  1C1L-.NT  (  ZERO  ORDER) 

C.  FN  -  SIGNAL  F  Oil'  I  f  K  Cni  FF  I  C  i  F  NT  S  (VECTOR) 

r.  FAVf.O  «•  AVER  A  (it  l>  SIGNAl  C.OEFP  (ZERO  ORDER) 

f.  FAVG  *  AV(RMrt>  SIGNAL  FOURIER  COEFF  (VECTOR) 

C  RHO  =  NUT  AT  1 ' ‘ S  RADIUS 

C  SIC,  =  STANOA-O  DEVIATION  uF  P.LUR  CIRCLE 

G  WI  s  WIDTH  (F  DETECTOR  ( X- 0 i REG T I  ON) 

C  h  e  HEIGHT  (it-  0:T  EOT  ufc  (  Y-  0  IP  EC  T  ION  ) 

G  Al  *■  IGVIRSC-  C/>K  LENGTH  OF  rt  AC  K,  GKOUN  D  C X-L1RDCT I  ON) 

G  Rr  =  INVFkSF  CD<--  UM.i'H  CTF  t,  AC  KGR  CUNT  (  Y-D I  KfcC  T  1  ON) 

G  Pfl  c  S 1  AN  IT  AF  I.'.  DEV  1ST  i  0  N  ji  POINT  JOG  ERROR 

C  XO  =  DrTFCTC-,  COORD  I NA  T  f  a-DIRLCTION 

G  YO  =  OF  TFT.  TOR  C  TORN  1  N  A  Tr  Y-DIRICTIUN 

G  kX  *  POINT  TARGET  COORDINATE  X-DIRCCTION 

G  RY  =  POINT  TARGET  COORDINATE  Y-DI  RECT  1  ON 

G  R .  k  AN  Cl.  F  *  PfJl  AG  COORDINATES  OF  POINT  TARGET 

G  NS  -  NilMIlFR  Of-  T  I  RST  CCfFFIClrNT  DESIRED 

G  NUM  =  NUM'TES  OF  LAST  GC.EFf  ICIFM  DESIRED 

G  IAVE-1  "SIGNAL"  WILL  CALCULATE  AVERAGE  SIGNAL  COEFF 

f.  I  A  Vf  =0  "SICNAl"  WILL  CALCULATE  POINT  TARGET  COEFF 

G 
C 

im  ICIT  kFAtvn  i  A  —  ri  *  D  -  /  ) 

DIMENSION  A (  F  A (  12  01 

CCMPl  EX'  16  F N (  1^01 , FSO, FAVG( 1J0)  . FA VGO 

COMMON  /DT  PR  k HO  .  S  1  G,  a‘  1,  H .  AL  ,  oE  «  RX  » RY  »  XO  »  YO 

COMMON  /A  VPR  M /  ■UIM  ,  1  AVt  ,?0 

COMMON  /START/  NS 

COMMON  /  ANG1  /  kANGL 

REAO  ?.RHl,SIG,r,' I  ,H,Al  ,6c 

R  F AO  ? , XO. YO. R, ANGL 

READ  3 . NS .NUM , I  A VE ,PG 

P.  ANGl  -  Q  .0  1  7ADOO  '  ANGL 

RX*R-OCC.S  (  .ANGl  > 

R  Y  =  R '  DS I N( RANG  LI 

RANGL  *  DAT  AN?(  k  /  - YO.RX-XO) 

C At  L  SUMGK 

CALL  fl  K. GF  IN  (  A  1.  F  A  •  A  A  J 

C  Al  L  S  IGNAL<  FN.FNO) 

WRITE  (9)  RhO.S IC. w l ♦ h, AL . 3E.RX. RY. XO. YO. PQ 

I  A  Vl=*l 

IF(  lAVE.EO.l)  A  A  =  0 . 00  0 
IFUAVF.FO.i)  RY=O.ODQ 
GA(  t  SUMC.K 

GAll  S IGNAl(FAvG.FAVGn) 

WRITF  (V)  F  MO ,  F  N  ,  A  A  ,  A  rA 
WkITF  (9)  FA  VGO.  F  A  VG 
1  FORMAT  (  10  IS  1 

FORMAT! FD1 0.0) 

FORMAT <  3  1 S,  ID  10.0) 
format  lar  ioo ) 

STOP 
END 


lie 


nooo  on  noon  noon  noon  oonooooonoooooooo 


SUBROUTINE  6KGR  N ( V  P , V  1 


SUBROUTINE  BKGRN  COMPUTES  THE  COEFFICIENTS  OF  TH' 
BACKGROUND  CORRELATION  FUNCTION  UP  TC  120  ORDERS. 
HIGHER  ORDERS  MAY  EASILY  QL  COMPUTED  BY  ONLY  CHANGING 
THE  DIMENSION  STATEMENTS. 


ZP  »  VALUES  OF  RADII  FOR  DIFFERENT  ANNJL I 
2  *  WORKING  VECTOR  FOR  RADII 
MAXR  =  MAXIMUM  RADIUS  OF  OUTER  ANNULUS 
NI21  *  DIMENSION  0^  IP 

VP  *  WORKING  VECTOR  IN  FIRST  PART  AND  RETURN  VFCTOR 
IN  THE  SECOND  PART 
V  «  ZERO  ORDER  COEFFICIENT 


IMPLICIT  RGAL-9  (A-HfO-2) 

DIMENSION  VP (  1 )  *  ON ( SO  *  1 20 ) 

DIMENSION  DI(50),Zd(20),?.  (50) 

0 1  MS-  NS  1  ON  R  (  2  A  )  ,  W  (  2  4  )  .  T  H<  2  0  ) 

COMMON  /  A  V  P  R  N;  /  NUM  ,  I  A  Vt  ,  PfJ 

COMMON  /  DTP- M/  R  HO ,  S I  G,  W  1 ,  H ,  AL  ,  b E  ,  R  X,  R  Y ,  X 0 ,  YO 
COMMON  /BKPR  M/  A , B  ,  C  »  D 
COMMON  /NF-ICST/  N  F  I  R  S  T 
DATA  MAXK/4/,NR/l/,NIZl/20/ 

DATA  2 '5/0. 00  0,1. 2  DO,  2*  4  00 .5. 5  DO,  8. 6  00,1  1  .800,15.00,20. 
IDO, 25. 000, 30. 600, 40. 00 0.50. 000, 60. 000, 70. 0D0, 80, 0D0, 

2 90.000, 10 0.0  DO, 15 0.0  DO,  200  .0  00, 4  00.000/ 

NFIRST=0 

PI  *  0 ARCOS (-1.00) 


CALCULATE  CONSTANTS  FCR  INTEGRA! 

CONST  =  AL<BE FSIG‘»2/PI/Pl 
A=W I / 2 . 00/ S  IG 
B  =H/2 . D  0/ S  1G 
C  «=  (  AL-"  SIG)' R? 

D  «  !Bc*.‘S  JG)**2 
E  *  RH0/S1G 


SET  ANNULI  FOR  SPECIFIC  PARAMETERS 
DG  ISO  I  =  1 , N I Z  1 

150  zm-zpm/E 

1  «  0 

IF( Z(NI Zll .GT .MAXR >GD  TC  210 
DZ  *  PI /NR/E 
MMM*  50- NI Z1 
DO  200  1=1, MMM 

2  <  I  *N1  Z  1  >  =  2  «  N I  Z 1  )  ♦  I*  OZ 
1F(?( I ♦ N I Z 1 ) .GT.MAXRIGD  TO  210 

200  CONTINUE 


MAX  =  NUM5EF.  OF  ANNUL! 

210  M AX= 1 +N I  21-1 


PRINT  4 ,RHO,SIG,  H, Kl , AL , BE , A , B , C  ,  0 ,  E , DZ 
PRINT  4  .  ( 2 ( I  1,1=1, MAX  ) 

DO  300  1=1, MAX 


COMPUTE  RADI  I, WEIGHTS  AND  ANGLES  FOR  EACH  ANNULUS 


1 1  7 


c 

c 


CALL  WAA(£(I),2(  I ♦  1 ) » P  t TH » W ,NPQ ) 


NF I RST  =  1 
D1 ( 1 )  =  0.00 
00  lit?  NN-l.NUM 
115  DNU.NM)  «  0,  DO 


C  COMPUTE  FUNCTION  FOR  EACH  RADIUS 

c 

00  260  K«-l,NPO 

CALL  3F  SSL  ( K  ( <  )  ■'  E ,  VP  ,  VJN  ,  MUM) 

RV=VJN-  >  2  /  0  F  X  P I  F  (K  )*‘2  ) 

EP=0l:XPI«<K)  +  *2> 

DO  120  NN=1,NUM 
120  VP(NN)=VP(N.N)*#2/=P 

KECcO. DO 
C 

c 

C  COMPUTE  FUNCTION  FOR  EACH  ANGLE 

C 

00  270  J=1,NPQ 

27  0  RcCeRr.C*r-UUMK)  <  DCOS(TH(  J)  )  ,  A I  K )  *  0$  1 N  I T  H  (  J  ))  1 
V- 0 .DO 

DO  130  NN=1,NUM 

13  0  ONI  I  ,NN»  -DN(  1  ,NN )  +  REC<  VP(  NN) W I  K > 

280  DItll  =  01(1)  ♦  RtCVRV^W(K) 

300  CONTI  N  J  F 
C 

c 

C  SUM  value  FROM  FACH  ANNULUS 

c 

DO  400  T  =  1 . M A  X 
J  =  MAX  +1-1 
400  V  =  V  +  01 (J) 

V=V*C3NST 

PRINT  3 ,N,V,DI(MCX  1 
N  601  =  0 

PRINT  13.NBQ.V 

PRINT  4 , ( DI (  1  ) , I =1 , MAX  1 

00  500  NP=*1,NUM 

VQMeO.OO 

VPINP) -0. DO 

00  60  0  1=1  ,MA  X 

•1  =  M  AX  +  1-1 

VPINP) "VP ( NP) +DNI J,NP) 

1FIDNI J . NP ) . GT, VGM)  VOM-CNI J,NP 1 

600  continue 

VPINP)  =  VP  (  NP  1  ♦'CONS  T 
PRINT  13, NP, VPINP) 

PRINT  1  4  ,  NP,  V  QMf  DM  MAX ,  NP  1 
500  CONTINUE 

£  CALCULATE  SUM  FOR  INTEGRATION  CHECK 

C 

SUM*V 

DO  700  J « 1 , NUM 
I  =  N  UM  +  1  -  J 

SUM?  SUM  +  2 . 00s  VP  I  I  I 
700  CONTINUE 

PRINT  12,SUV 
RETURN 

1  FORMAT (101 5) 

2  FORMA',  (  of  10.0) 

3  FORMAT ( 21 10, 1P2F 15.6) 

4  FCRM/TI  lP8r:15.6) 

12  FORMAT ( 5X, G20. 1 3  1 

13  FORMAT  <5X,  14 ,6X ,020.13) 

14  FOKMAT( 5X,I4,6X,G20.13 ,6X,G20. 13) 

16  FORMAT  (  ;■>(  1P020.  12)  ) 

END 
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FUNCTION  FU(X, Y  ) 


FUNCTION  FU  15  AN  AUXILIARY  R'JUY  INC  ^UR  DKC.RN  THAI 
COMPUTES  THAT  PART  CF  THE  IN1EGRANU  NUT  RADIALLY 
SYMMETRIC. 


IMPLICIT  RCAL’-G  (A-H.C-ZI 
COMMON  /  UK  PRM/  A,B,C.,D 

P 1  *  1 . DO 

I  F  <  X.NE  .  0.  DO)  P)=DSIN(A-'X)/(A»X) 
1F(Y.NC  .O.DOI  P 1  *  P  1 :  I'S  1  N  (  EK:  Y  )  /  (  8M  Y  ) 
FU»PlH'2/(C<-Xf.: 2)  /(D«-Y*<2) 

R  ETURN 
END 
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SUBROUTINE  SIGNAL! REC ,RECO) 


SUBROUTINE  SIGNAL  COMPUTES  THE  FOURIER  COEFFICIENTS  OF 
A  POINT  TARGET  DR  AVERAGED  COEFFICIENTS  OVER  A 
GAUSSIAN  POINTING  ERROR 


ZP  =  V  A l U  Fc  OF  RADII  FOR  DIFFERENT  ANNULI 

Z  =  WORKING  VECTOR  FOR  RADII 

MAXR  =  MAXIMUM  RADIUS  CF  OUTER  ANNULUS 

N I Z 1  =  DIMENSION  L-F  ZP 

RFC  =  RETURN  VECTOR  OF  COEFFICIENTS 

RECO  -  ZERO  ORDER  COEFFICIENT 


IMPLICIT  RE  &  L  8  <A-H,0-Z> 

Cl  MENS  ION  BN (120) 

01  ME  NS  I  ON  01(50), DNJ ( 50 , 1 2 0  I , ON ( 50 , 120 ) , Z < 50 ) , Z P  (  20 } 

01  MENS- 1  ON  R  (  2  4  )  ,  W(  24)  ,  THl  24) 

COMPLEX" 16  RFC! 1 ) , VFC( 120 ) , RECC, VS* F,  WP 
COMMON  /A  V  c'  R  M  /  NUM,IAVt,PO 

COMMON  /DTPRM/  RHG , SI G, Wl , H, AL , BE ,RX,RY , XO, YO 

COMMON  / S I P  R  M /  A , B , C , PX , QX 

CCM.MCN  /ME  I  R  $T/  NE1RST 

COMMON  /START/  NS 

COMMON  /ANGL/  RANGL 

DATA  ZR/ 0.00  0,1.  2D0, 2. 4 DO, 5.500,6.600,1  1.000,15.00,20. 
IDO ,25 .000, 30.000, 40.000, 50.000, 60. 000, 70. ODO, 80. ODD, 
290. ODO,  100. 000, L5 0.0 DO, 200. ODO, 400.0  DO/ 

DAT  A  M  a  XF  /  7  i  ,  \R  /  1  /  ,  N  I  Z  1  /20/ 

NFIPST  ~0 

VSsDCKPLX! O.CO.-l .CO) 

FI  -  OAPCOSt- l.CO) 

P 1 2  =  P 1/2  .t j 
KKK-0 


CALCULATE  CONSTANTS  FGR  INTEGRAL 
GAMA=S I G 

I  F(  1AVF  .EC.  1  )  GA,MA=C$QRT<  SI  G-  T2  +  PC**  2) 
SG  =  DS  CRT  (  2  .  DO  )  *'  GAMA 
A  =  W I  /  S  G 
B  =  H  /  S  G 

CCNST  =  A'<B/PI/PI 
C  =  1.00 
PX=“XO 
OX  =-Y0 

IF (IAVE.EQ.O)  P  X  =  P  X  +RX 
IF(  IAVT.EO.G)  5X*QX+=Y 
PX=PX<  DSORT  (2  .CO  )  /GA'.' A 
QX=QX*’DS0RT(2  .  CO)  /  GAMA 
E  =  2.D0-PH3/SG 


SET  ANNULI  FOR  SPECIFIC  PARAMETERS 

OZ  =  0  .  DO 

DO  150  I=1,NIZ1 

Z  (  I  )=  ZP ( I ) / ( 2.D0"E ) 

I  =  0 

IF(ZCNIZl)  .GT. MA  X  R ) GO  TO  210 

C2=  15  .  DO”  PI/E 

M M M = 5 0 -  M  ll 

Dr'  200  I=1,MMM 

ZU+NI71)  =  ZINIZl  )  +  3 A  D  Z 

If  (  Z  i  I  +  N I  Z 1 ) . G  r . MA  XR ) GO  TO  210 

CC'NT  I  NO  5 

PRINT  1  7  ,  Z  II  «-N  I  Z  1  ) 


noon  onon  noon  oon 


MAX  =  NUMBER  OF  ANNULI 


210  MAX=I +NI Zl-l 
GO  TO  212 

211  MAX* I 

212  PC  1  NT  4  ,RH?  ,  SIS  ,Wl  ,H,  AL  ,BE  ,  RX.RY  .XO.YG,  PX,  QX 
215  P  F<  1  NT  4  »  (  Z  (  1  ),  1  =  1,  MAX) 

220  DO  300  1=1, FAX 


COMPUTE  RAD! 1 .WEIGHTS  AND  ANGLES  FOR  EACH  ANNULUS 

CALL  WAA(Z(  !),?.<  I  +  1J,R,TH,W,NPQ1 
NF I RS  T  =  1 
NPR=NP0 
OKI)  =  0 .  DO 
DO  115  NN=NS .  NUM 
D\J<  I  ,NN )  =0.00 
115  Cf  (  I ,  NN  )  *  0.  DO 

1 F ( KKK. 10.1)  GO  TO  300 
DC  230  K-l.NPQ 
RK2  =  R ( K  i  C  +  ? 

IMRK2.GT.174)  GO  TO  <85 
EP=0LXP (RK2) 

IF(FP,GT.1.D  50)  GO  TO  285 
CALL  BE  SSL  (P  (K ) . E  »BN , VJN ,NUM) 

R  V  =  v  J  N / FP 
PECG=F (  Rl  K) » T  H . 0 ) 


COMPUTE  FUNCTION  FOR  EACH  RADIUS 

RFCC-RECQ*K\M  W  (  K  ) 

PP*RECO 

D1  (  I ) *C 1 ( I )i PP 

DO  130  NN  =  N S .  NJM 

BN ( NN)  =BN(NN)/fo 

REC  (NN)  =F(  p  U) .Th.NN) 

REC(NN)=REC(NN)<BN<NN)>rw(K) 

PP=REC (  NN* 

CN(  1  .NN  )  =  ‘ > N (  1  .NM  +  PP 
PP- REC ( ON  )  VS 
ON  j  (  1  ,NM“LLJ(1  .  NN  )  +  °P 
130  CONTINUE 
280  CONTINUE 
300  CONTINUE 


SUM  VALUE  FROM  EACH  ANNULUS 

VZ  =  0  .  DO 
VI  =  O.CO 
V2  =  O.DO 
DO  400  1=1, MAX 
J  =  MAX  +1-1 

IF (DABS (DI( J) ) .ST. LABS! VZ 1 )  VZ=DI ( J) 
1F(DI ( J ) .lE.O.OO)  VL  =  Vl+01 ( J ) 

400  1  F  (  D I  (  J  )  .  G  T  O.DO)  V2  =  V2+DKJ) 

Rfc  CO- CCMPl X( vl  +  V2, 0.00) 

Rf  C0=  c  ECC->  CONST 
N6Q  =  Q 

PRINT  4, ( DI (  1  ) , 1  =  l.MAX  ) 

PRIM  1  9  ,  N  P  C  ,  F  L:  C  0 
PRINT  14.NBD.VZ.DI (MAX) 

DO  500  NP=NS i NJM 
VOM  -  O.DO 
VI  =  O.DO 
VI  J  =  O.DO 
V  2  =  O.CO 
V2J  =  0.D0 
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00  600  1*1, MAX 
J=MAX* 1  -I 

1F10NI  J,NP).l  T.0.00)  Vl=Vl+DM 
I  F(DNJt  J,NP)  .IT  .0.00)  VlJ=VlJ+PNJ(  J,NP) 

1 F(DN(  J ,NP) .GE. 0. DO)  V2=V?+DN(J  ,  N  P ) 

IF(CNJ( J,NR) .GE.0.C0)  V?J=V?JfDNJ( J,NP) 

I F(DAbS (DN( J, NP)  ) . GT . DABS ( V Q.M ) )  V QM= DN( J  »  NP I 
600  CONTINUE 

PEC(NP)*DCMPLX(V1+V2,  V1J«-V2J) 

RECUP)  =  F  E  C  ( NP) *  CONST 
PRINT  1 9,N°,REC( NP) 

PRINT  14, NP,VQM, DN(MAX,NP ) 

500  CONTINUE 
725  CONTINUE 

VS  =  OCMP  LX ( 0 . DO,  1 >  DO ) 

DO  800  NP* NS , NUM 
800  PEC(NP)=FEC(,\|P)*VS+‘‘NP 
I F ( IAVL.EQ.0 )  GG  TO  960 
SUM  =  R  ECO 

I F ( NS  4 \ E . 1 )  SUM*  0 . ODO 
PRINT  1 9 , N60,  RECO 
C 

c 

C  CALCULATE  SUM  FOR  INTEGRATION  CHECK 

C 

DO  900  I  =NS ,  NUM 
PQ=REC(  I  ) 

PRINT  19tI.RECm 
900  SUM=SUM+2.D0*PQ 

PRINT  18,  SUM , SUM J 
RET  UPvN 

960  CONTINUE 

V  S  =  DC HP  L  X ( 0. DO,-  1  .DO) 

PRINT  19,NPQ,RECC 
SUM=RECO 

IFtMS.NE.i)  SUM*  0 . DO 
SUM J-0 . DO 
DO  950  NPsNStNUM 

VFC(NP)=P.EC(NP)  vDCMPLX(DCOS(NP-RANGL  ) »  0  S I  N C  NP*RANG L  »  ) 
PRINT  19, NP, VEC(NP) 

PP= VEC ( NP  ) 

SUM=SUM  +  2  .DO*PP 
PP=VEC ( NP) * VS 
SUMJ=  SUMJ  +PP 
950  CONTINJE 

PRINT  1 8 , SUM , SUM  J 
RETURN 
285  KKK= 1 

GO  TD  300 

1  FORMAT (10 15) 

2  f-C'RMAT (  8E  10.0) 

3  FORMAT (2110,  IP  25  15.6) 

A  FORMAT ( 1P8EI5. 5) 

5  FORMA  1(215) 

6  FORMAT ( IP2E15. 6) 

12  FORMAT ( 5X,G20.13) 

13  FORMAT (  5X, I  4, 6X, G20. 1 3) 

14  FORMAT  <5X , 14 ,6X, G20.12.fcX, G20. 13) 

L 5  FORMAT (  2X, I  3 ,2X, 1 PE15.o ,2X, 1PE15.6,2X, I  3, 2X.1PE15 .6, 2X 
1 ,  1  PE  1  5 . 6  ) 

lfc  FORMAT  (31 1 R 0 2 0 . 1 2  )  ) 

17  FORMA  T(  5X,  »m;.x  \,jt  LARGE  ENOUGH  Z*  Ax  =  ’  ,  c  1 5 . 3  ) 

18  FCRMAl ( ?X , • SUM* • , G20. 13, • S JMj* 5 , G20. 13) 

19  FORMAT (5X,I4,6X,G25.13,G20.13) 

END 
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FUNCTION  F  (R.TH.NN) 


FUNCTION  F  IS  AN  AUXILIARY  ROUTIN'.  C0R  SIGNAL  THAT 
COMPUTES  THAT  PART  OF  THE  INTEGRAND  NOT  RADIALLY 
SYMMETRIC 


IMPLICIT  REALMS  (A-H.O-ZI 
COMPLEX'  16  F 
DIMENSION  T H <  1  I 

01  ME NS I  ON  AAE (24 I , LBE (2  A ) , A AO (24  >  ,  DBO( 24  I 
COMMON  /S1PRN/  A,B,C,PX,QX 
CATA  NPP/24/ 

J.)e(-1»YVNN 
AA=  0.  000 
BB=0 .000 

IF(NN.GT.O)  GO  TO  5 
DO  1  J  =  1  *  N  P  0 
X  =  RADCOS(TH<  J  )) 

Y=R>'DSI  N(  Tti(  J  ) ) 

DSX=OS  INI  PX  X  ) 

0SY=0S1  M  QX  Y) 

DCX=OC  >  S(  PX"  X  > 

CCY  =  DCOS ( QX '  Y  I 
P)  =  1.U0 

I  FlX.NF  .0.00)  P  UPS  INI  A*  X)  /'(  A«  XI 
IF(Y.NF.O.DO)  PI  =  P1-:  DS  INI  BTY  )/  (  BT  Y  ) 

A  AE  (  J  I  -  OC  X  ■  OC.  Y»  P  1 
BBE(  JUOSX  "OSY'-Pl 
A  API  J)  =  -PCX  DS  Y  '  PI 
B8G(  J  U  ~P  S  X'.  DC  Y •' P  L 
1  CONTINUE 

5  HI  J.LLT.OI  GO  TO  10 

00  20  J  =  1 f  NP  9 

AA=AA+AAEl  Jl  OCQS  INV  THIJ  )  I 
6S=BS+BEE(  J)  '••OS l  NlNNrTHUH 
20  CONTINUE 

F®OC'vPLX(  AA,  B6) 

IF(NN.cQ.O)  F =DCMPLX( AA,C. DO) 

RETURN 

10  DO  40  J=1 ,NPR 

A  A=  AA  +  ft  AC  <  Jl  *  0  SIM  NM  TH  (  J  5) 

83=6F  +  o EC.  I  J  )’ GCCS INN  THt J  )  ) 

40  CCINT1 NUF 

F  =  PCMPL  XI AA, 66  > 

RETURN 

END 
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SUBROUTINE  WAA{  R  l ,  R  2  ,  P  ,  TH  ,  W  ,N) 


f. 

C 

f.  SUBROUTINE  WAA  COMPUTES  THE  RADIAL  PO I  NT  S  ,  WE  IGHTS  AND 

C  ANGIES  AT  WHICH  THE  INTEGRAND  MUST  BE  EVALUATED, 

f. 

r. 

C  N P  «  DEGREE  OF  GAUSSIAN  LA 

f.  R  «  RADIAL  POINTS  TO  RE  JSED 

C  TH  *  ANGLES  TO  BE  USED 

C  We  WEIGHTS  TO  BE  USED 

C  NFIRST  PREVENTS  CCMPJT ING  G-L  POINTS  MORE  THAN  ONCc 

C. 

c 

IMPLICIT  REALMS ( A-H.O-A) 

DIMENSION  R(il,W(ll,TH(ll 
RFAL  *8  RC(  24)  ,,-JCt  2-4) 

COMMON  /NFIRST/  NFIRST 
DATA  NP/24/ 

P 1 cOARC  OS (  —  1 .DO) 

PINePI /?. DC/NP 

IFINFIR ST.NF.OiGO  TO  99 

CALL  GL024 (RC.mC.Q .00. 1 .DO ) 

DO  10  1 =1 „NP 
10  TH<  1)M  I-.500)  •  P  IN 
NFIRST  =  1 

99  CONTINUE 
N  =  NP 

T  =  R2  -  Rl 
TW  =  TR  <R2  *-  Rl  I 
F  =  TW.PI/NP 
DC  100  1  =  1. NP 
RX  =  Rl 
R  X  *  RXSRl 
RY  =  n? 

R  Y  a  R  Y=rR 2 

R  (  I  )  =  DSORTCftX  ♦  RCIUMRY  -  RX)) 

W II)  =  WC ( I ) L  F 

100  CONTINUE 
RETURN 
END 
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noOr>o 


St.'RKOUT  I  Nc  Gt.024  (X,  A.C.  D) 


GAUSS)  AN-LFGENDER  INTEGRATION  ROUTINE  (24  POINTS) 


Or.URLfc  PRECISION  DMC.OPC 

DOUBLE  PRECISION  C.D.Xl  24). A<  24), XXI  12»,  AA(  12) 
DATA  XXm,AAU).XX<2).AAm.XX(3),AA(3),XX{4),AA(4),X 

1  X  X ( 6  J  «  A  A  1 6 )  . X  X ( 7 )  ,AA<  7)  .XX  (8)  .  AA(tt)  ,XX  (9  )  ,  AA(9  ).XX  UO 

2  xxim.AAi  m.xx(  l?).aai  12)/ 

* .9951 872199970210601 7999740970-0, .12341229799987199540 
*.474 , 26555971 309*48 19 8 391 9 93  -0 . . 2 853 1 3 8862893366 3 181 

938274 55 200273 2 7 6 fa 52 3 6*9001 <u-0, .44277*38617419806163 
4 . 886*15  52  7  004401  03421 315434 140-0,  . 5429 6 5849  1 34 3o7 8 0 74o 
* .82000 1983 97 6902 42 1953949 872o0-0.. 73 3464 8 1 4 11 080 30 57 3* 
< .7401 24 1 VI 57855* 30*2438 28 10 j 11- 0. .86190161531953275917 
*. 64809  3  651 5364  7  5  5692 5? 49 5  78 o9 0-0  .. 976 1 3 6 5 2 1 04 11 36 8 82 69 
*-  .  5*542  l  *71  38883*5  356  38275  61  720-0, .  107*44  27011596563  4  78 
<’  .433  79  3  507626  0*5  I364fc7084c  3150-0.  .11550566805372560135 
*■.3150*2  679646163  37  43  66  79329 1  3  0-0  .  .  1 2  1  o  7  04  72  92  75  03  39 1 20 
* . 19  111  8  86747  361  6  309 13  8639  8  2  070- 0. . 1 2 5b 3 7 *  56346 82 8 2 96 1 2 
4. 6*056 6 9? 86260562 6 08 50* 3 08 2 62 D- 1 , .1279381953467521569/ 
OUC  =  .  5U0MD-C) 

DPC  =  .  500;*  (  0«-C  I 
00  2  1=1.12 
N 1  *  25  -  I 

XU)  «=  -DMC'XXI  I  )  ♦  OPC 
X<  N1 )  =  CMC-  XXU  )  -t-  DPC 
A  ( I  )  =  OMC-r-tAI  I  ) 

2  A  INI)  =  use. AA(1) 

RETURN 

END 


m 


SUBROUTINE  Bfc.SSL  (  AG.RF  ,BZ  ,NUM> 

c 

c 

C  Su-sRnuTlNF  fir  S  S  L  CCMPU7FS  BESScL  FUNCTIONS  CF  THE 

f.  RECEIVED  ARGUMENT 

C 

f. 

f.  AG  =•  ARGUMENT  OF  RISSFL  FUNCTION 

f.  BF  *  KFlUP.si  VcLTOF.  OF  dFSScL  FUNCTIONS 

f.  BZ  *  ?  F  K* )  OFOEK  Hi  SSEI  FUNCTION 

C  NUN  -  HIGHEST  CmGfV  OF  BESSEL  FUNCTION  TO  bE  RETURNED 

f  NY .NX  =  OROfRS  CALCULATED  BY  ASYMPTOTIC  EXPANSION 

C. 

c 

IMPtlOIT  REALMS  (A-H.fl-2) 

DIMENSION  RF3( I0OO) ,6F ( 1) 

N7= ) 000 
KTsO 

IFIAG.I  T„0)  K  I  s  1 

IHAG.l  T.O.DO)  AG*  DAB  S  (  AG ) 

X=AG 

NY*)F 1 XI SNGl  ( AG)  ) 

1 F(NY .1 E. 1 )  NFs lb 
IFINY.F-O.  ?)  N  F  ~?  0 
iriNY.FQ. 3)NF  =  2<. 

]  F  (  NY  .  G  F  .4  .  AND.  NY  . LT  .6  )  N  r  *  Y 
I  F(  NY.GF.  6.  AND.  .mY.  LT.  I  0)  NF  NY 
1FINY.GF.I0.AnD.NY.lT.1S)  NF  =  A*  NY 
1  FINY.GF.IS.ANO.NY  .LT  .29)  NF=3»-NY 
1 FINY.GF .29.AND. NY.LT. 8G»  N  F  =  2  -  N  Y 

I  F<  NY  .GF.  60.  AN.).  NY  .LT  .  ISO)  N  F  =  )  f  I  X(  1 .  S<  F  L  0  A  T I  NY  )  ) 

I  F  (  NY.GF  .1  SO.  AN'J.NY.LT  ,33  0  )  NF  =  1  F  1X1  1  .<.*  FLUAl  <  NY  )  ) 
IFINY.GE.330)  N  F  =  I ? I  X (  1.23  FLOAT ( NY) ) 

lf(lrlXIFLOAT(NFI/2.).LT.IFIX(FL&AT(NF+l)/2.))  \F=NF+1 
IF(NY.LT.I)  NF  =  3 
N  Y  s  N  F 
NXsNY-1 
C 

C  SET  DROcRS  >  NY  TC  ZERO 

C 

DC  1  I=NY.W 
B  F  3  (  I  )  =  0.  no 
l  continuf 

c 
c 

C  CflMPUTF  ASYMPTOTIC  EXPANSION 

c 

DO  2  MU  =N  X  » NY 
Z*X/DFL OAT (mu ) 

APU=DFlC AT (MU) 

A  »  OSUkT(1,C3-Z*»2) 

A  l  si  .SOu-.  I  01  0G(  I  1  .DOt-A)  /  l  )-A) 

DEL  =  DFXP(?.D0/3.C0#DLCG(Al)) 

,v.  1  =  OF  XP  I  2.  DO/ 3.00-0  LOG!  A. MU)  tD  LOG  (DEL)  ) 

CALL  TABLE  (XI , A  IV.AI PV ) 

01  =  GEXP(  1. 00/3. OO-OLOG  (AMU)) 

02  =  CcXP( 5.00/ 3.00’Dl OGlAvj) i 

BG  s  -S  .00/(40.  f'O’  OF  L  *  *2  )  ••(  L  .DD/DSCRT  (  DEL  )  )  -^  (  S-DO/(  2h. 
100''  A  *  *  3  )  -  1 .  Oi)  0  /  <  8.000-  A  t  ) 

FI  =  OtX° ( 0.2 500  DlDG< 4 .00* DEL I -0 . 500 >D LOG  ( A) ) 

RF3(MU)  =  F I *  ( At  V/01 ♦ A  IPV  FB0/D2 ) 

?  CONTINUE 
NY?  =NY-  2 
C. 

r. 

C  USF  RFCURSlVc  EOUATIDN  to  generate  orders  <  NX 

c 

no  i,  jo.ny? 

N  =  ,-.t  2-1+1 

B  F  3 { \ )  =  (2. DO- D FLOAT  (N  +  l  )/X)  -  o  F  3  (  N  ♦  1  )  -  3  F  3  (  N+  2  ) 


^26 


f>psfSo  r»r>r»o 


4 


DCn  ORDER  BESSEL  FUNCTIONS  ARE  ODD 
IFtKT.FG.l)  BF(  I)*BF(Ii*l-l.DO)«--»  l 
77  CONT Ijlp^0j04fiF3(l)/X-BF3(2J  )  *  SNO  RM 
RETURN 
END 


3  2  7 


SURF  OUT  I  Me  TABLE  {  X,  A I  V.  A I  PV) 

C 

c 

C  SURkfiUT  INF  TABLE  IS  AN  AUXILIARY  ROUTINE  FOR  BESSL 

C.  THAT  COMPUTES  AIRY  FUNCTIONS 

C 

f. 

IMPLICIT  R F A L » 6  (A-H.O-Z) 

DIrtFNS  ION  A  I  ('30  >.  A  1  P  (  SO  ) .  Art  I(  15  )  ,  AOIPC  15) 

DATA  A1  /4  +  0.0  00,  .  562?  80  DO.  4*0  .0  00  ,  .  560*6200,  V-0 .000* 
1.5567241/0,4'  u.OOO,  .55705800,4*  0.  ODO , .  55  54  5600 , 4-  0 .  ODD, 
2 .5535 W  DG.4' 0 .OUO, .55242100, 4>0 .00 0.  .55096000,  4+0.  ODO, 
3.54458400,4;  0.000,  .  54823000/ 

OAT  A  A  I  P/4.  0.000,  .  5o 6 6  730 0.  4,  0.  ODO,  .  56944  80  0,4*0,000, 
l. 57192  7  00.46  0. 0 00. .57432 000,4 <  0 .000, . 57663500 , 4 *0 . 000, 
2 .  5766  7800.  4-0. 00  0.  .  58  1  05600 . 4^-  0 .  ODD ,.  583  1  74  00 ,4*0. 000, 
3. 5852.16  lO.  4-  O.uOO,  .  5872-, 500/ 

DATA  Ah  1  /  4  •'  ( i .  0  0  0  .  .5482  30  00.  .5456  36  00,  .54  316000, 

1 . 54064', 00.  .63861  800.  .  53648900..  53444800,  .  53248800, 
2.530601 00. .52678300, .52702700/ 

DATA  Art  1  p/46  0. 000 ..  587245  00  .. 59  1  12000, .5  94  823  00, 
i  .  648  372  DO. .60  178  20  0,  .  60506  600. .  6  082  3500, .61130500 , 
2.6142  75  DO. .61715600, .61995400/ 

L  =0 

Y  *  l)FXP(DLnr,<  1 .5001-1. 5D0»DLDG<  X)  ) 

W  =  DEXP(DinG(2.D0/3.D0)*l .5DO*DIOG(X) 1 
L  =  1 

1  f  (Y.GT  .1  .500)  GO  TO  14 
IFtY.GT.O, 44995 900)  GO  TO  1 
L  =  2 

W  1  =  10 .00  *  Y 
I W 1 = 1 F  I  X ( SNGL (Will 
W?=Wl+  1  . 

I  W5-  1  F  I  XI  SNGt  (  w?  )  ) 

I  Z  M  W 1  +  1  0 

IFdWl.LT.IW2)  IZ*IZ+5 
I Y= 17-6 

I  F(  I  Y.  FQ.O)  C.G  1  0  13 
Gfl  TO  a 

1  Y  =  1  0 . 00  «  Y 

I Yt I F I X<  SNGL ( Y 1 ) 

I F (  IY-Y  1  2.4,6 

2  I  Z  =  I  Y+  l 

GO  TO 8 

4  A  I V  =  A  1  (  lY  ) 

A  I P V=  A  I  P(  I  Y) 

GO  ro  1  5 
6  I  7*  lY 

I  Y= I Y-  l 

8  0P=  Y-FL  CAT( I Y ) 

I F ( L. EO  .2 )  0P=  OP/ 5 .00 
IF(L.Fo.l)  GO  TO  12 
A  1 V  =  (  1 .0  0- DP  )  ‘  A  l<  I  Y  )+DP>A  I  (  I  Z  ) 

A  I  P  V  s  (1. DO-OP)  «AIP{  IY)+0P*AIP(  IZ) 

GO  TO  19 

1?  AW  =  <  1  .DO-OP)  vABII  1Y)+0P*ABI  d  Z) 

A I  PV  =-  U.OO-DP)  *A3IP(  lYi+OP”ABIPdZ) 

GO  TO  1 9 

13  OP  *  Y / 6.00 

A  1  V  =  I  1  .  DO-OP)  0.56415003+0° -A  I(  IZ  ) 

A  I P V  =  ( l . 00-DP ) +  0. 564 15000  + Op^AIPi  IZ) 

19  AW  =  OtXPI-Dl  jG12.D0)-0.25D0^0LCG(X)-W  +  0LGG(  AIV)  ) 

A  I PV  =  02xP(0.?500J  OLOGCX)  -* +Dl OG ( A  I  Pv J  ) 

A  I P V  =  A  IP  V  < - . 5D0) 

RFTIJ3N 

14  P  3  1 T  30 

30  FOs-iAT  (  5X.  *H^V3  EXCEEDED  TABLE  VALUES') 

RCTIKN 

ENO 
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r>oon  r>r»r»r*  r»nr>r»  or>nnr>r»r>r'>r>r>r> 


SUHROUT IN?  SUMCK 


SURROUT  INF  SUMCX  COMPUTES  T HF  ANALYTICAL  EXPRESSION 
FOR  THE  INFINITE  SUM  CF  COEFFICIENTS 

SI  SUM  =  INFINITE  SUM  OF  POINT  TARGET  fIAVE-1)  OK 

AVFRAGFO  TARGET  HAVE-0)  COEFFICIENTS 

BKSUM  v  INFINITE  SUM  uF  BACKGROUND  COEFFICIENTS 


IMPLICIT  REAL*  8  (A-H.C.-7.  J 

COMMON  /DTPRM/  RhG . SI G , Wl ,H , AL , BE . RX ,RY , XO. YO 
COMMON  / A  V  P  K  M /  NUM.IAVE.PO 
COMMON  / ANGl /  SANGL 
W  =  WI 

CONST  *0 .750000 
GAMA-DSORT  (SI  G»  F2  +  PO*-*-2) 

1HIAVF.F0.0)  GAMA  =  SI  G 
DIV-2 .OO.DSURT ( 2. 00) -GAMA 
I F ( lAVE.EO.OJ  GO  TO  10 


COMPUTE  SUM  CF  AVERAGED  COEFFICIENTS 

A-DERF  (  (  2..C0  M  XO+RUO  H  W)/01  V) 
fl - OFF  F (  (2  .1)0  <•  (  X  0  ♦  RHO )  —  W  )  /  0  1  V  ) 

C  =r-ERF  (  <  2  .  00  OH  YO)  +H)  /DI  V) 

0  =  f>FRF<  (2.000H  YO)-H)  /  D I  V  ) 

4  SI  SUM-CONNl  -  (  A-b)  <• ■  (C-D) 

R  1=R<  W.  AI  .  SIG) 

R2^R(h, PE. SI G) 


COMPUTE  SUM  OF  BACKGROUND  COEFFICIENTS 

RKSUM*4.00*AlME*Ri*R2/(*n>H)H2 

PRINT  5 . PKSUM, S I  SUM 

RFTUFN 

TO  CONTINUE 


COMPUTE  SUM  GF  POINT  TARGET  COEFFICIENTS 

A-GER  F  (  r?,Ono*  (  XO  +  RHO--DCOS  4  RANol  l-RX  )  +  W  )  /  C I  V) 
R  - 0 F R  F  <  (  ?  .  00  0*  (  XOf  RnC*  DCCS  (  SANGL  )  -RX  )  -«  )  /  C  I  V  ) 
ODFRF(  (  ?  .OOM  YO  +  fiHO*  OS  I  N(  R  ANGL  )  -A  Y)  +»-i)  /C  l  V) 

0  -OF  RF  (  (2  .  DO*  <r  D  +  hhG'  OS  INI  F  ANGL  )-RY  )-rl)  /  DI  V  ) 
GO  TO  4 

5  FORMAT  <2X,'RkSJM='.015,8,10X.'SI SUM-* ,D 15.6) 
END 
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1 


r»r>r>r»r> 


FUNCTION  R  ( II  «Z  «  S 1  G I 

FUNCTION  R  IS  AN  AUXILIARY  ROUTINE  FOR  BKSUM 


lMPllf.1T  R£M*8  (A-H.o-Z) 
p  1  eUAkCOSI-1 .DO) 

fI«-<(U/(P.OO*SIG)  )•*•■*?) 

E1«DKXP(F1) 

•  F ? *  < Z S  IGI  'w? 

E?*Df  XPIF2) 

F  3)  =  7  f  U 
F3*DFXP  (Fill 
FAtDF  XPI  -  c31  ) 

Xl  =  l)/  l  ?  .DO’  S  1G) 

Y 1  =1  .DO-OFPHX1  ) 

X?«=7*  SIG+U/IP.DO^SIG) 

Y?*1.00-DFRF(X?I  r 
X3*Z*  S  I  G-U/  (  2  .  CO"  S  IG) 

Y3“  1. DO-DEAF  (  X3) 

XA  =  Z+S  IG 

R»UJ/?ZRv!?Wfl.D0-Yl-<SI3*2.D0/(DSQRT(Pl  )*U))* 
R*M(E?/(?  .D01>Z<  *3 J  I  >  <  E3AY2  +  EA<-Y3-2.D0*YA) 
RFTURN 
END 


(1 . 00-E 
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